Pion production in nonequilibrium Chiral Perturbation Theory 
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We apply the formalism of Chiral Perturbation Theory out of thermal equilibrium to describe explo- 
sive production of pions via the parametric resonance mechanism. To lowest order the lagrangian 
is that of the Nonlinear Sigma Model where the pion decay constant becomes a time-dependent 
function. This model allows for a consistent nonequilibrium formulation within the framework of 
the closed time path method, where one-loop effects can be systematically accounted for and renor- 
malized. We work in the narrow resonance regime where there is only one resonant band. The pion 
distribution function is peaked around the resonant band where the number of pions grow exponen- 
tially in time. The present approach is limited to remain below the back-reaction time, although 
it accounts for nearly all the pion production during the typical plasma lifetime. Our results agree 
with the analysis performed in the 0(4) model. The space and time components /J 1 (t) are also 
analyzed. To one loop /° 7^ unlike the equilibrium case and their final central values are lower 
than the initial ones. This effect can be interpreted in terms of a reheating of the plasma. 

PACS numbers: ll.10.Wx, 12.39.Fe, 25.75.-q, 11.30. Rd 
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I. INTRODUCTION 

Nonequilibrium Quantum Field Theory has attracted 
a considerable attention over the past decade, partially 
motivated by the experiments on Relativistic Heavy Ion 
Collisions (RHIC). The ultimate goal of such experiments 
is to describe the properties of the QCD phase diagram 
and the Quark-Gluon Plasma (QGP) from the observa- 
tion of the final hadronic spectra. The Relativistic Heavy 
Ion Collider at BNL, which is already running, will reach 
energy densities high enough to confirm the very promis- 
ing results obtained at CERN (SPS) during the last few 
years, indicating the existence of the QGP [GJ. For a re- 
cent update of both theoretical and experimental results 
in this area see 0. A lot of theoretical effort has been 
put in trying to understand the various properties of the 
QCD phase diagram in different ranges of temperatures 
and densities in thermal equilibrium. In addition, there 
are many aspects of the nonequilibrium behaviour of the 
plasma which are not fully understood. In the standard 
picture of the collision |3|,[ij the plasma formed in the 
central rapidity region cools down very rapidly, reaching 
approximate local thermal equilibrium. During the sub- 
sequent expansion, the temperature scales of chiral phase 
transition and deconfinement are crossed and hadrons are 
produced. In this regime, observables depend only on 
proper time approximately, in the central region. The 
expansion goes on until the final freeze-out of hadrons. 
The typical plasma lifetime during which nonequilibrium 
effects are important and most of the final hadrons are 
produced is about 10 fm/c. 

A possible scenario to explain the observed final 
hadron distributions is that where strong fluctuations of 
the pion field are formed during the chiral phase transi- 
tion, giving rise to the so called Disoriented Chiral Con- 
densates (DCC) M. These were suggested originally as 



misaligned vacuum regions where the chiral field is point- 
ing out in a different direction in isospin space from that 
where the vacuum expectation value of the pion field van- 
ishes. Ideally, one would have misaligned regions of ob- 
servable size, with only neutral (or only charged) pions. 
If such regions were formed, one would observe large clus- 
ters of pions emitted coherently from the plasma as the 
pion field relaxes to the normal vacuum [jfj. This kind 
of behaviour is indeed observed in Centauro events in 
cosmic ray experiments 0. A clear signal for DCC for- 
mation has not yet been observed in RHIC experiments 
|| . However, it seems that one has to measure higher 
order pion correlation functions in order to identify a 
pure DCC signal which is not masked by other effects 
P,|l0[ so that the search will continue at BNL. Other ob- 
servable consequences of DCC-like configurations would 
be an enhancement of dilepton and photon production 
JTlIjOI and a modification of the effective 7r° — > 77 ver- 
tex |13|. In any case, pion production within the energy 
scales of the chiral symmetry provides a natural frame- 
work for hadronization |L4j . 

Thus, one should be able to describe nonequilibrium 
phenomena such as large pion production, from the mi- 
croscopic theory governing the relevant degrees of free- 
dom. At the energy scales where the chiral symmetry 
plays a predominant role (below 1 GeV) QCD is non- 
perturbative and one has to use an effective lagrangian 
which describes satisfactorily the microscopic meson dy- 
namics. Such a theory must incorporate the QCD sym- 
metries and the chiral spontaneous symmetry breaking 
(SSB) pattern. In this picture the Nambu-Goldstone 
bosons (NGB) are the lightest mesons (it, K, rj) and the 
masses of the light quarks are meant to be treated pertur- 
batively. One possible choice is simply the O(N) model 
where the fundamental fields are N — 1 pions and the a, 
and the potential has the typical SSB shape. However, 
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one should bear in mind that this model becomes non- 
perturbative in the coupling constant at low energies so 
that it is imperative to perform alternative expansions 
such as large N. On the other hand, the O(N) model 
shares the QCD chiral symmetry breaking pattern only 
for N = 4, so that it is not able to incorporate kaons and 
etas. 

An alternative approach is an effective theory built as 
an infinite sum of terms with increasing number of deriva- 
tives, only in terms of the NGB fields. The Nonlinear 
Sigma Model (NLSM) is the lowest order action one can 
write down in this expansion. Higher order corrections 
come both from NGB loops and higher order lagrangians 
and can be renormalizcd order by order in energies, yield- 
ing finite predictions for the meson observables. The 
unknown coefficients, which encode all the information 
on the underlying theory, absorb the loop infinities and 
their finite part can be fitted to experiment. This frame- 
work constitutes the so called chiral perturbation theory 
(ChPT) Jl^,[l6| which provides a well-defined perturba- 
tive expansion in terms of p/A x , where p stands gener- 
ically for any meson energy scale of the theory (masses, 
external momenta, temperature and so on) and the chi- 
ral scale A x ~ 1.2 GeV (see 0,ff§ for a review). One 
of the many advantages of this scheme is that it can be 
extended from the SU(2) chiral symmetry (only n fields) 
to SU(3) with K and rj. The ChPT formalism has also 
been applied in thermal equilibrium to analyse various 
properties of the low temperature meson gas [p~9[— ^T|| . 

In the context of nonequilibrium chiral dynamics, two 
possible scenarios for pion production and DCC forma- 
tion have been proposed: the first one takes place in the 
early stages of the plasma evolution. Roughly speaking, 
after a very rapid cooling the chiral field is at the top 
of the classical potential in the chirally broken phase. 
As the field rolls down, long wavelength modes grow ex- 
ponentially (spinodal instabilities) and this behaviour is 
responsible for the enhancement of DCC's. There have 
been several approaches in the literature to implement 
this idea in the 0(4) model |6 p^j23[| . Typically, the pion 
distribution function is peaked at low momenta, being 
different from a thermal distribution |23j whereas the 
pion densities and DCC sizes predicted are around ~ 
0.2 fm -3 and 1.5-2 fm respectively. The second sugges- 
tion is based on the parametric resonance mechanism 
[ pj and inherits the idea from inflationary reheating p4| . 
The analysis in the spinodal regime shows that the time 
it takes for the field to roll down to the bottom is very 
short compared with the total plasma lifetime. Thus, 
in the parametric resonance approach, the a field is os- 
cillating around the minimum of the potential in a later 
stage of the plasma evolution. Those oscillations transfer 
energy to the pion modes, giving rise to pion solutions 
exponentially growing in time via parametric resonance. 
Typically, the unstable modes develop in bands in mo- 
mentum space and the more important resonance band 
is centered at k ~ m a /2 P,p5|. The DCC sizes in this 
approach can be as large as 5 fm ]2(| and recent cal- 



culations show that strong charge-neutral correlations in 
parametric resonance can be used to identify a pure DCC 
signal Q. Furthermore, the reheating process yields pre- 
dictions for the final hadronization temperature compati- 
ble with the observations |l4|] . One must stress that both 
approaches are complementary and in fact the initial con- 
ditions needed for the parametric resonance correspond 
to the final stage of the rolling down solution. A very 
detailed analysis of both regimes in the context of the 
O(N) model can be found in [f27| . 

The purpose of this work is to explore pion production 
in parametric resonance within the ChPT framework, as 
a complementary analysis to the 0(4) model. We will 
show that within this formalism one can also describe re- 
gions where the number of pions and the pion correlator 
grow exponentially. The main advantages of the ChPT 
approach are that one can follow a consistent perturba- 
tive treatment which is renormalizable order by order and 
that it can be extended to three flavours. Besides, one 
is dealing only with NGB fields, although we will show 
how the pion production can be understood in terms of 
the a field evolution in the 0(4) model. This method is 
best suited for the stage of the plasma expansion where 
the system is well into the broken phase of the chiral 
symmetry. This is precisely the regime where parametric 
resonance takes place. 

We will build on a previous work pq] , where we have 
analysed the extension of ChPT out of thermal equilib- 
rium. In that work it has been shown that the power 
counting and renormalization program can be consis- 
tently implemented also at nonequilibrium. In turn, the 
present analysis will provide a particular example where 
it will be shown explicitly how the chiral power counting 
and renormalization program work, yielding predictions 
for physical observables. The key idea is to make use 
of the derivative expansion consistently implemented in 
ChPT in order to study the system not far from equilib- 
rium. For that purpose, the nonequilibrium dynamics is 
encoded effectively in the parameters of the model. To 
leading order and assuming a spatially homogeneous sys- 
tem, we let the pion decay constant be time dependent. 
This function acts as an external force on the pion degrees 
of freedom. It is important to bear in mind that a self- 
consistent treatment should amount to incorporate the 
full hydrodynamics of both the fluid and pion modes p9[ . 
In the present approach we will concentrate only on the 
influence of the expansion on the meson dynamics. This 
is a similar situation as considering a QFT in the pres- 
ence of an external curved background space-time |30"|| . 
In that case, it makes sense under certain conditions to 
ignore the back-reaction effect of the matter fields in the 
metric. Similarly, we will see that for the time scales rel- 
evant for pion production it is reasonable to ignore those 
effects and treat the influence of the expanding plasma 
as external. 

The paper is organized as follows: In section || we will 
review the nonequilibrium ChPT and its relationship to 
curved space-time QFT, which will be crucial in what 
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follows. The parametric resonance ap pro ach in its sim- 
pler version will be discussed in section III , while sections 
IV and [y] will be devoted to analyse the effects of para- 
metric resonance in two different observables: the pion 
decay functions and the pion number respectively. The 
latter is the most relevant observable as far as pion pro- 
duction is concerned whereas the former will allow us to 
estimate the time scale when the back-reaction effects 
become important as well as the final temperature by 
that time. In both cases we will calculate up to one-loop 
in ChPT, paying particular attention to renormalization. 
Besides, in section fvj the definition of particle number we 
will use and its relationship with the energy-momentum 
tensor are discussed. Our conclusions are summarized 
in section [v|. We have included two appendices collect- 
ing some useful results about the Mathieu equation and 
curved space-time QFT respectively. 



II. NONEQUILIBRIUM CHIRAL 
PERTURBATION THEORY 

A. The NLSM out of thermal equilibrium 

The chiral lagrangian to lowest order is the nonlinear 
sigma model (NLSM), which contains two derivatives of 
the fields. In the chiral limit, where the mass of the light 
quarks is set to zero, the NLSM only contains one energy 
parameter / . To lowest order, / is nothing but the pion 
decay constant / = f„ ~ 93 MeV (/ ^ f n to higher or- 
ders) . Consistently, our nonequilibrium model will be the 
NLSM where / is replaced by a time-dependent function 
f(t). This function acts as an external field encoding the 
time evolution of the system as, for instance, the expan- 
sion of the plasma formed after a RHIC in proper time. 
As for the initial conditions we will assume that the sys- 
tem is in thermal equilibrium before some initial time 
(t = for convenience) at the temperature T; = fij ■ 
This is an important simplification from the point of view 
of the nonequilibrium path-integral formulation. In fact, 
one can formulate the generating functional for the real- 
time Green functions by extending the time arguments 
to the contour C in the complex plane showed in Fig- 
ure 0. This is an extension of the Closed Time Path 
technique [ jsT) for nonequilibrium field theory, where the 
inclusion of the imaginary-time leg is a consequence of 
the choice of equilibrium initial conditions p^ . Thus, 
f(t < 0) = /, and the function f(t) represents an exter- 
nal force switched on at t — and driving the system 
out of equilibrium. Note that we choose that departure 
to be instantaneous and then f(t) cannot be analytical 
at t = 0. Nevertheless, it is natural to expect that these 
restrictions on the initial conditions do not influence very 
much the physical results for longer times, especially if 
the system is not far from equilibrium. 

Thus, our starting point will be the following nonequi- 
librium NLSM action: 
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FIG. 1. The contour C in complex time t. The lines C\ 
and C2 run between ti + ie and tf +ie and tf — ie and ti — ie 
respectively, with e — > + . 



S 2 [U] = [ d 4 x Lj9. tr d^(x,t)d^U{x,t) (1) 
Jc 4 

where J c d 4 x = J c dt J d 3 x. We will restrict here to the 
case of two light quark flavors (i.e, the NGB are only 
the pions) and hence U{x) S SU(2). Besides, we will 
be interested only in the chiral limit i.e, massless pi- 
ons. Thus, we are not including any explicit symmetry- 
breaking term in the action. Note that the action (|l|) 
is chiral invariant (U — ► LUW where L, R are constant 
SU(2) matrices) by construction, which will play an im- 
portant role in what follows. As it is customary, U(x) is 
parametrized in terms of the pion fields 7r a (a?, t) as: 

U{x,t) = -^-{ [f(t) - TT 2 (X, t)] I + iT tt 7T B (2, t)} 



where I is the 2 x 



2 identity matrix, r a are the Pauli 
matrices, n A = Tr a ir a and TT a (ti — ifli) = n a (ti) is the 
equilibrium (KMS) boundary condition, with ti < 0. As 
we will show below, the physics does not depend on the 
choice of t i . 

It is useful to recall that the NLSM can be viewed to 
lowest order as an 0(4) model subject to the constraint 
a 2 + 7r 2 = f 2 . The same applies to the nonequilibrium 
case and therefore we can think of f(t) as the lowest order 
expectation value of the a field in the 0(4) model, which 
has not reached its equilibrium value yet. 

With the NLSM one can predict the very low energy 
behaviour of pion observables such as the pion decay con- 
stant or the pion scattering amplitude, which agree with 
the current algebra predictions. However, to go beyond 
the lowest order, one has to consider pion loops. Thanks 
to Weinberg power counting theorem [I5j we know that 
the loop diagrams are of the same order as lagrangians 
with more derivatives of the pion fields. In fact, it can be 
shown explicitly that the (undetermined) coefficients of 
such lagrangians absorb all possible UV divergences com- 
ing from the loops and hence one gets finite and scale- 
independent predictions for the pion observables. 

In the present case, we need an extra ingredient to 
the power counting, namely, the time derivatives of the 
function f(t). We will take 
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and so on. In this sense, one remains close to equilib- 
rium. The rest of the power counting is the usual one, 
i.e, every derivative of the pion field is 0(p) and every 
pion loop introduces an extra 0{p 2 ). Other than being 
subject to the conditions (Q), we will let f{t) be arbi- 
trary. However, in this work we shall discuss how f(t) 
can be chosen consistently with physical results such as 
pion production. 

Let us now expand the NLSM action to lowest order in 
the pion fields. Using the boundary conditions, we can 
integrate by parts and write, to second order in the pion 
fields, 



Sa[n]=~ [ d 4 xn a (x,t) [u + m 2 {t)} n a (x,t) 
2 Jc 



(3) 



where m 2 (t) = -f{t)/f{t). That is, the model accommo- 
dates a time-dependent pion mass term, without break- 
ing explicitly the chiral symmetry (unlike a physical pion 
mass term). This is an important feature of this model, 
since it suggests an interesting connection with a QFT 
in the presence of an external curved space-time back- 
ground. We will discuss this point in detail in section 
[IC. Before that, let us analyse the main properties 



of the leading order two-point function (the propagator) 
arising from the above NLSM action. 



B. The leading order propagator 

The two-pion correlation function is iG ab (x,x') = 
(TcTT a (x, t)ir b (x', t')), where Tc indicates time ordering 
along the contour G. Note that, from isospin invariance, 
we can write just G ab (x,x') = 8 ab G(x,x'). Besides, by 
spatial translation invariance the two-point function de- 
pends only on x ~ x' . This is not true for the time co- 
ordinates due to the nonequilibrium time dependence in 
the lagrangian. 

There are four different types of propagator depending 
on the relative position of t and t', namely, G 11 , G 12 
and so on [Q. Instead of writing all the combinations 
explicitly, we shall keep the condensed notation of time 
ordering with Tc defining the natural extensions 9c, 8c 
and so on. Thus, from (||), the leading order Go(x,x') is 
a solution of the differential equation 

{d 2 + k 2 + m 2 {t)} G (t, t', k) = ~5 c {t - (4) 

where Go(t,t',k) is the Fourier transform in the space 
coordinates only and k 2 — \k\ 2 . 

As for the boundary conditions, thermal equilibrium 
for t < means that we have to impose KMS boundary 
conditions at the imaginary-time leg in Figure [j] [Q . That 
is, defining 



G(t,t',k)=G > {t,t',k)9 c (t-t') 

+ G < (t,t',k)e c {t' -t), 



(5) 



the KMS boundary conditions read 

Gq (ti — t ,k) — G^(ti, t , k) 



(6) 



The general solution to the differential equation (Q) 
with the boundary conditions (||) can be constructed for 
all the branches of the contour in terms of two particular 
solutions hi t2 (t,k) to the homogeneous equation [[33): 

{df + k 2 +m 2 (t)} h % (t,k) = i = 1,2, (7) 

such that their Wronskian 

W(t, k) = hx(t, k)h 2 (t, k) - hi(t, k)h 2 (t, k)^0 (8) 

It is important to remark that the general solution 
Go(t,t' , k) must be continuous and differentiable in the 
time coordinates so that it is uniquely defined. Thus, 
in our case we demand that hi(t,k) for t > and their 
first derivative match the equilibrium solutions at t = 0. 
Since m 2 (t < 0) = 0, two independent equilibrium solu- 
tions are given for t < by 



h?(t,k) 



-ikt 



'2k 



h e 2 q (t,k) 



1 



2k 



(9) 



which we have normalized so that W eq (t,k) = 1. It is 
not difficult to see that continuity and differentiability at 
t = imply that W(t, k) = 1 also for t > 0. On the other 
hand, since m 2 (t) is real, if hi(t) is a solution to ([?]), so 
it is h 2 (t) — ih*(t), where W(t) = 1 if it matches the 
equilibrium solution at t — 0. 

We shall be dealing here with real time evolution for 
positive time coordinates and therefore, unless other- 
wise stated, we will be interested in G 11 only. In that 
case, we will suppress the "11" superscript for simplicity. 
Nonetheless, it should be borne in mind that in the loop 
integrals there are contributions from all the branches of 
the contour pS}] . 

Thus, the solution for Gj 1 is given by |33]| 



hi{t,k)hl{t',k)6{t-t') 
h\{t, k)hi{t' \k)9(t' -t) 
n B {k) [/n(t,fc)ftl(t',fe)4 



hl(t,k)hi(l/,k)] 
(10) 



where t and t' are both positive and the boundary con- 
ditions at t = imply 



ht(t = 0+,k) 



2k 



h 1 {t = Q + ,k) = x j'- (11) 



where the dot means d/dt. The dependence with the ini- 
tial temperature appears through the Bose-Einstein dis- 
tribution function 



n B (k) 



1 



oftfe _ I 



(12) 



In different parts of this work we will need the two- 
point function evaluated at the same space-time points: 
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G (t) = G (x,x) 



(27T)° 



r(^)( 47 r)^ Jo 



dkk d - 2 G (t,t,k) (13) 



where the equal-time correlator in momentum space 
reads, from (jig): 



iG (t,t,k) = [l + 2n B (fc)] |/n(*>*)| S 



(14) 



Because of the loss of time translation invariance, Go(t) 
is a time-dependent quantity. Besides, it may be UV di- 
vergent and, therefore, we will use dimensional regular- 
ization (DR) with d the space-time dimension, which is a 
suitable regularization scheme as far as chiral lagrangians 
arc concerned |lq [111. 



C. Renormalization and curved space-time 

Once we have defined our noncquilibrium power count- 
ing, we can apply ChPT to calculate the time evolution 
of the observables. In doing so, we must pay special 
attention to renormalization. The fact that there is a 
time-dependent mass term indicates that there can be 
new time-dependent infinities in the chiral loops. For in- 
stance, in standard ChPT with a nonzero pion mass, the 
tadpoles renormalizing the pion propagator to lowest or- 
der yield the usual infinities proportional to in DR. 
These infinities are absorbed by two counterterms pro- 
portional to m 2 in the fourth order lagrangian p6| . We 
expect similar divergent contributions here proportional 
to m 2 (t). However, we are working in the chiral limit 
and therefore we are not allowed to introduce the above 
mentioned counterterms. Otherwise we would break ex- 
plicitly the chiral symmetry. Hence, we should be able 
to construct the most general fourth order action, which 
in particular has to include new terms (and hence new 
low-energy constants) to cancel those extra divergences, 
preserving exactly the chiral symmetry. 

There is a natural way to find this 0(p 4 ) lagrangian, 
using a very fruitful analogy: the action ([I]) can be writ- 
ten as a NLSM on a curved space-time background cor- 
responding to a spatially flat Robertson- Walker (RW) 
metric, with scale factor a(t) = f(t)/f. For that purpose 
it is more convenient to work in terms of rescalcd fields 
n a (x,t) = n a {x,t)/a(t). Hence, U =I(/ 2 - n 2 ) 1/2 I + 
iT a n a )/ f and we can write the action (Q) as 




ggi MV tid ll U^{S,t)d v U(x,t) (15) 



S 2 [U] 



where the metric g^ v is nothing but the spatially flat 
RW metric with line element ds 2 = a 2 (t)[dt 2 — dx 2 ] in 
conformal time t |5(j| and g = det g = — a s (t). 

With our chiral power counting, it is straightforward 
to assign the chiral order of the covariant tensors con- 
structed from the metric. For instance, the Ricci tensor 



R^ v is 0(p 2 ) and so is the Ricci scalar R — g^ v R pv , and 
so on. Explicit expressions for these tensors and other 
useful results for this metric are collected in Appendix 
[b|. An important point in this formulation is that we are 
considering the so called minimal coupling of the mat- 
ter fields with the metric. That is, we are discarding 
possible couplings between the pion fields and R{x) to 
0(p 2 ), such as £R(x)tr(U + f/t). The reason is that we 
want to preserve chiral invariance, which would be bro- 
ken by those terms |34[]. Thus, in this language, m 2 (t) in 
(^) represents the minimal coupling with the RW metric 
preserving chiral invariance. 

Therefore, we have a systematic way to construct the 
nonequilibrium lagrangian to any order. We just have 
to include all possible terms consistent with the chiral 
symmetry, contracting indices covariantly with the met- 
ric g^ iv (x). In particular, to 0(p 4 ) it reads 

S i [U,g,I{\ = [ d 4 x^[C 4 [U l9 } 



- (LuRg^ + LiaR^ixdpUidvU] (16) 

where £4 [U, g] stands for the standard (equilibrium) la- 
grangian ]l6| ] with indices raised and lowered with the 
g^" metric and the rest are new 0(p 4 ) invariant cou- 
plings with R(x) and R^ v (x) in the chiral limit. These 
are the new terms we need, where in and L\i are the 
new coupling constants. We are following the notation of 
Ipif , where this lagrangian was first considered to stud; 
the energy-momentum tensor of QCD at low energies 
The same lagrangian has been used in the context of pion 
hard exclusive production |56| . In J3J] it has been shown 
that in order to cancel the one-loop infinities, Ln has to 
be renormalized as 



ljJL° 



-4 r 



6 16tt 2 



11/ 

— --(l0g47T-7+l) 



(17) 



with d the space-time dimension, 7 the Euler constant 
and n the renormalization scale. L{\(/i) is finite and 
depends on fi so that the combination in the r.h.s of 
( |l7| ) remains scale-independent. On the other hand, 
L12 is finite. Their numerical values can be obtained 
from the experimental information on the QCD energy- 
momentum form factors. They yield L12 ~ —2.7 x 10~ 3 
and L r n (fi = 1 GeV) ~ 1.4 x 1CT 3 . We will use these 



*The terms with Lu and L\% in ( |lq ) differ in a global 
sign from those in H|. The reason is that we are following 
here the convention in ^] for the Riemann tensor, namely 
Rp^s ~ — • • •; where r«L are the Christoffel symbols 

(see Appendix |b|) whereas in jS) the convention for Rp^s i s 
reversed in sign. For the same reason, every term propor- 
tional to Ln or L12 here has its sign changed with respect to 
those in H. 
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same values here since their possible nonequilibrium cor- 
rections are of higher order in our analysis. 



Thus in our case we have to replace in (16) the RW 
metric. After expanding in the 7r a fields and partial in- 
tegration we have: 

Sa[tt, g] =-\f d 4 OT° [h{t)d 2 t - f 2 {t)V 2 + m\{t)] ir a 



+ 0(tt 4 ) 



(18) 



with 



hit) 
h{t) 
m\(t) 



■12 



{2Ln + L\2 
(6Ln + L 12 ) 



, f(t) 
' f 3 (t) 

/(*) 
Pit) 



[/(*)]* 
1 Pit) 

f(t)f 

nt) 



h(t)Rt) + h{t)f(t) i .. 

W) + 2 /lW 



(19) 



The above lagrangian should take care of the nonequi- 
librium infinities we might find in the pion two-point 
function. As far as this work is concerned, these are 
the only infinities we will have to renormalize. 

In the following sections we will concentrate on a par- 
ticular case for f(t) (or the scale factor if we use the 
curved space time terminology) which is of physical rel- 
evance as the simplest approximation producing a large 
number of correlated pions. In addition, this example 
will allow us to test explicitly the cancellation of the 
(new) nonequilibrium infinities appearing in the observ- 
ables considered. 



III. PARAMETRIC RESONANCE AND PION 
PRODUCTION 

A. The parametric resonance approach 

It is clear that our approach will be useful in a stage 
of the plasma evolution when the departure from equi- 
librium is of the same order as the meson energies. 
Hence, as far as pion production is concerned, we are 
in the parametric resonance regime. Let us briefly re- 
view some of the ideas behind parametric resonance in 
the 0(4) model In the last stage of the field 

evolution, the a field is oscillating near the true vac- 
uum and those oscillations have relatively small ampli- 
tude [^3|. Following a semiclassical approach, one can 
split the a field as a(x,t) = ao(t) + 5a(x,t) where <7o(t) 
is a time- dependent homogeneous classical background, 
solution of the equations of motion to leading order in 
the amplitude, whereas 5a includes next to leading or- 
der corrections and quantum fluctuations. One can pro- 
ceed perturbatively around the classical solution ao(t). 
In a first approximation, quantum fluctuations of both 
the a and the pions can be neglected. Thus, if the a 



field is oscillating around the potential minimum a = f 
and the amplitude of the oscillations q is small, one 
can solve the equations of motion perturbatively in q. 
To leading order the equations of motion for the a and 
the pions decouple from each other and one simply gets 
<7o(t)/f = 1 — (9/2)[sin(rn CT t + if) — simp] where m a is 
the a mass, ip is an arbitrary phase and we have chosen 
the initial conditions so that the field is at the bottom of 
the potential for t = 0, consistently with our choice of the 
initial equilibrium state. Neglecting the pion correlations 
is equivalent to state that (tt 2 ) <C f 2 ■ If the first order 
solution is inserted into the equations of motion to next 
to leading order, the pions satisfy a Mathieu equation. 
The importance of this equation is that it has solutions 
exponentially growing in time for certain bands in mo- 
mentum space. This is the essence of the parametric 
resonance mechanism, which even in this simple classical 
picture is consistent with hadronization [ fl4"| . The para- 
metric resonance idea is directly imported from reheating 
and preheating in inflationary cosmology p4| ] , where the 
small q approach is called the narrow resonance limit for 
reasons to become clear below. 

The approach described in the previous paragraph is 
the crudest one can follow in this context, although it re- 
produces the main features of parametric resonance. One 
can refine it in several ways. First, neglecting the pion 
correlations but keeping the NLO terms in the a ampli- 
tude leads to a Lame equation instead of the Mathieu 
equation p7| , p6[ . As it is emphasized in |2^J, the differ- 
ence is not only quantitative, but the resonance struc- 
ture is also different. As for the importance of pion cor- 
relations, one must bear in mind that this is a matter 
of time scales. Such correlations grow exponentially in 
time until eventually they reach the same order of mag- 
nitude as the a term. As it is customary, we will refer 
to that time scale as the back-reaction time, so that, 
typically (7r 2 )(£##) ~ f 2 . The name is again inherited 
from Cosmology where the back-reaction describes the 
modifications of the metric or the inflaton field due to 
quantum fluctuations of the matter fields |24|,|3C]]. It is 
clear that for t > tsR,, pion correlations must be in- 
cluded self-consistently in the dynamics of the a field. 
In O(N) models this usually requires numerical simula- 
tions. For instance, in the large N limit it has been shown 
that, when the back-reaction is properly accounted for, 
pion correlations are responsible for the damping of the a 
field from tsR onwards and, more importantly, this dissi- 
pation stops the exponential growth in the pion number 
p7| , j37| . A different story though is that one can interpret 
that process as thermalization, or in other words, that 
the final particle spectrum is thermal. In fact, that is 
not the case when pion amplification occurs p3| . Thus, 
in practice, all the interesting physics associated with 
pion production in parametric resonance takes place be- 
fore the back-reaction time. Finally, a word must be said 
about the quantum corrections of the a field, which can 
be treated scmiclassically |J. In the narrow resonance 
approach, the evolution of the pion fluctuations is influ- 
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enced only by terms linear in <7o(t) and the dynamics 
of the cr fluctuations is not important for pion produc- 
tion, since the width of the resonance band for the a is 
negligible with respect to that of the pions. Thus, the 
essence of the exponential growing of pion fluctuations is 
not changed qualitatively by including quantum correc- 
tions in 5a. 

At this point, let us establish the connection with our 
present ChPT approach. Our philosophy will be to work 
out consistently the simplest choice for f(t) yielding para- 
metric resonance. For that purpose, it is useful to com- 
pare with the 0(4) model. In the oscillatory regime (with 
small oscillations) one starts by keeping only the lead- 
ing order a ~ oo(t). In that limit and for small oscilla- 
tions the NLSM is nothing but the 0(4) model subject to 
7T 2 + a 2 = (Xq (i) . In other words, we should take simply 



/(*) = / 



1 - -sinMi 
2 



(*>0) 



(20) 



where we have chosen ip = so that 

m 2 (t) = -(qM 2 /2) smMt + 0{q 2 ). 

Note that with this choice, not only /(0 + ) = / but also 
m 2 (0 + ) = 0, i.e, both f(t) and m 2 (t) match their equi- 
librium values at t = 0. This w ill play an important role 
in the analysis in section III B . 

Thus, we would use the NLSM model (@) with f(t) 
in (20) if we were interested in describing pions out of 
equilibrium classically, for times where the back reaction 
is not important and the plasma is in the broken phase. 
Then, we can calculate pion observables using ChPT, 
where the pion fluctuations can be treated quantum me- 
chanically in a consistent fashion. As a matter of fact, 
there is no need to invoke the 0(4) model in the first 
place, since the NLSM is the lowest order action com- 
patible with all the symmetries, driven out of equilibrium 
through the time dependence in f(t). For that reason, we 
have replaced the mass of the sigma by M, an arbitrary 
mass parameter. We expect it to be in a range compat- 
ible with m a but the advantage of our present approach 
is that we do not have to worry about the uncertainties 
related to the a particle []. We will discuss below what 
range of values for M is reasonable in order to reproduce 
the right order of magnitude in pion production. 

Let us be more precise now about the smallness of 
q in our approach. According to our previous discus- 
sion about the noncquilibrium chiral power counting, we 
should demand at least that qM 2 = 0(p 2 ) and so on. 
In this way, all the 0(p 4 ) corrections will remain under 



^The a mass is not very well determined. It ranges between 
400-1200 MeV according to the latest PDG data Q. It is not 
even clear that one can describe it as a particle. For instance, 
it shows up in pion-pion scattering in ChPT as a rather broad 
resonance in the I = J = channel [tS9l. 



control, as we will see below. It is important to remark 
that in this work we shall restrict to one loop in ChPT. 
Going beyond that would require additional restrictions 
on the value of M. Before carrying on, we would like 
to summarize the assumptions and the limitations of the 
present approach: 

i) We are neglecting the possible back-reaction correc- 
tions to f(t) in j20|). This is valid for times below Ibr, 
when the pion correlations are of the same order as the 
leading order, that is, {~K 2 ){tBii) ~ f 2 - However, we will 
see in section IV that within our approach one can es- 



timate tsR by calculating the loop corrections to f w (t). 
Fortunately, as we said before, nearly all the relevant 
nonequilibrium pion production physics happens before 
that time. This limitation comes from the fact that we 
are treating f(t) as external, similarly to quantum field 
theory in an ex ternal curved space-time, as we have seen 
in section II C . A self-consistent approach, analogous to 



treat also the metric quantum-mechanically in a pertur- 
bative low-energy fashion would be very interesting but 
is out of the scope of this work. 

ii) We are assuming that the system is the late stage 
of the expansion, so that it makes sense to treat the am- 
plitude of the oscillations q as a small parameter (nar- 
row resonance approximation). That means we will only 
retain the leading order in q and thus ignore 0(q 2 ) cor- 
rections. This is consistent with ChPT to one-loop if 
qM 2 = 0(p 2 ). We will see that this simplification 
amounts, among other things, to consider the Mathicu 
equation for the pion modes. 

iii) We shall restrict to one loop in ChPT, in the chiral 
limit and for the SU (2) chiral symmetry. As commented 
before, the nature of the ChPT approach allows to ex- 
tend our calculations including quark masses and three 
flavours. 



B. Dimensional Regularization of the LO 
propagator. 

According to the previous discussion, the differential 
equation ([7|) becomes to leading order in q, the Mathieu 
equation: 



d 2 h\(z, k) 
ah 2 



+ (a(k) - 2qcos2z)/ii(z,fc) = (21) 



where z = Mt/2-n/4, a(k) = 4k 2 /M 2 and, without loss 
of generality, we will take q > 0. 

The solutions of the Mathieu equation are known and 
tabulated. We have collected in Appendix [A] some use- 
ful results about this equation. For our purposes, the 
most relevant feature is that it admits unstable solutions 
exponentially growing in time for certain values of the pa- 
rameter a. This is the simplest version of the parametric 
resonance mechanism. The instabilities arise in bands 
in k, centered at k n = nM/2, of width Afc„ = 0(q n ) 
(see Appendix [A]). Therefore, in the narrow resonance 
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FIG. 2. Profile of iMG (t,t,k) for T t = M/100 and 
q — 0.1. The instability band for this case lies roughly be- 
tween 0AM < k < 0.6M. 



use either the asymptotic or numerical solutions for the 
Mathieu equation. In the second, which is UV divergent, 
we replace the solution (|2^). It is clear that the piece 
proportional to the Bose-Einstein function is UV finite, 
since ns(fc) decreases exponentially for large k. In ad- 
dition, the integral of any power of k is identically zero 
in DR, so that dkk d ~~ 3 = — J^ a dkk d ~ 3 and we can 
absorb that contribution into the first piece in (^3|) [J]. 
Therefore, let us write for t > 0: 



iG (t) 



T? 



qM 2 




12 32tt 2 
1 



/(t)+iGT(*) 



iG (t,t, k)- — [l + 2n B (k)] 



, v sin Mi - %■ sin 2kt ] , 
o fc2 _^ 2/4 (24) 



approximation we will just neglect the width of all the 
bands but the first one. A typical unstable solution for 
the equal-time correlation function Go(t,t,k) has been 
plotted in Figure || around the first band for a particular 
choice of the parameters. The solutions typically oscil- 
late with an exponentially growing amplitude inside the 
unstable region. These will be the field configurations 
responsible for explosive pion production. 

Our next step will be to analyse the equal-time cor- 
relator Go (t) in parametric resonance, separating its UV 
divergent part in DR. Therefore, we need to know the be- 
haviour of Go(t, ij k) in ( O ) for large k. Since we are con- 
sidering only one resonant band in momentum space then 
Go(i,i, fc) is in the stable region for large enough k. For 
small q, the analytic solution in the stable zone is given 
in flA6|), Let us define k such that for k > fc (> M/2) 
one can simply take q — q/[2(a — 1)] in (A7). As ex- 
plained in Appendix |A|, it is not difficult to estimate nu- 
merically the value of kg for a given q. We have found 
that ko = M satisfies the above requirement within our 
approximation range for q < 1. Nevertheless, we have 
checked that our results do not depend on the choice of 
k as long as k > M. 

Thus, replacing the approximate solution ( A6) in (|A5|), 



solving for the coefficients A(k) and B(k) and replacing 
the result in(0) we find, to leading order in q, 



iG (t,t,k) = — [1 + 2n B (k)} {1 



qM 2 



4 k 2 ~M 2 /4 



M 

sin Mt sin 2kt 

2k 



(22) 



for k > ko. First, let us split the k integral in ( J13| ) as 

poo pk() poo 

/ = / + / (23) 

JO JO Jk 

The first piece is finite so that we can take d — 4 and 



where we have separated explicitly the equilibrium con- 
tribution which is just 



1 

2^2 



dkk- 



1 
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and the integral 



r/ n „, f°° „ sin2fci 

I(t) = M / dk— — 

K ' J ko k 2 - M74 



(25) 



(26) 



is finite for t > 0. It is important to point out that if 
we had chosen a different phase ip in /(t), as for instance 
f(t) = /[l — (g/2)(cosMi — 1)], we would have found a 
singular behaviour near t = 0. In fact, instead of I{t) 
above, we would have an integral which is finite for t > 
and d — > 4 but logarithmic divergent in the t — > + limit. 
In this sense, t acts as a natural regulator. This would 
not have been a limitation to our approach, since we are 
meant to observe the system for times such that the t = 
effects are unimportant. In fact, we have checked numer- 
ically that the influence of those terms is irrelevant for 
Mt > 1. The behaviour at i = is just a consequence 
of our non-analytic approach, where the nonequilibrium 
effects appear instantaneously and is a well-known prob- 
lem in nonequilibrium field theory. In fact, in p0[ it has 
been pointed out that it can be cured by a suitable choice 
of the initial state, which for a time-dependent mass term 
amounts to take m 2 (0 + ) equal to the initial mass. This 
is exactly what we have done with our choice of phase in 
(pp|), since m 2 (0 + ) = in the chiral limit. Our results 
confirm the analysis in pOfl from a completely different 
viewpoint, namely working in path integral within the 
ChPT framework in the DR scheme. 



'The above result can be understood formally in DR by tak- 
ing linwo JT fc d -7(fc 2 + m 2 ) = -k d ~ 2 /(d - 2). 
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The divergent part in ( pjj ) is given by 

sin Aft qM 2 f°° „ k d - 3 



r(^i)(47r)^ 



4 



dk- 



M 2 /4 



(27) 



We will proceed now to regularize this expression. We 
have H, 



Ud-3 



dfc 



TP 

2-fl 



k 2 -M 2 /4 4-d 

(fc 2 - M 2 /4) 



d 
2 : 



1,2 2 ; 4fc, 2 



d Aft 
o 



4-d 

which replaced in (|27|) yields: 
qM 2 



+ C)(d-4), (28) 



iGg*(t) 



32tt 2 



2(fc 2 - M 2 /4)— sin Aft 



i(log7r-7 + 2) 



+ 0(d-4) 



(29) 



As we will see in detail in section IV, we will be able 



to renormalize this divergence in the low-energy constant 
in so that the answer for the observables is finite and 
scale independent. This means that the regularization of 
the above UV divergence is consistent. 



IV. THE PION DECAY FUNCTIONS f n (t) TO 
ONE LOOP 

The pion decay constant acquires one-loop correc- 
tions in ChPT. Those corrections are finite and scale- 
independent once the contribution from the 0(p A ) la- 
grangian is taken into account |L6| . The same will happen 
in our nonequilibrium model, where the pion decay con- 
stant becomes a time-dependent function f n (t), which 
to tree level is just f(t). For us, the importance of the 
calculation of f v (t) is twofold: first, it will provide an ex- 
plicit check of consistency of our renormalization scheme. 
Second, it will help us to understand the time scales. In 
particular, the size of the loop correction will define the 
back-reaction time, when it equals the tree level contri- 
bution. For times well below that scale, our approach 
remains perfectly valid and yields predictions for observ- 
ables such as f-n{t) and the pion number. As a matter of 
fact, the same philosophy is followed at finite tempera- 
ture in equilibrium, where in the chiral limit fll9|]: 



■ JhT] - = J1{l -h 



(30) 



with fir ~ 93 MeV. Clearly, this result is valid only for 
temperatures below T* = v^ft- ^ 228 MeV. In fact, 
even though f-niT) is not a good order parameter for 
the chiral phase transition pl| , it should decrease as T 
approaches the critical temperature T c . Therefore, the 



one-loop result (|3C|) already reproduces the correct qual- 
itative behaviour and indeed it provides a reasonable es- 
timate T* for the critical temperature. The T 2 term in 
( |30| ) is nothing but the thermal pion correlator (tt 2 (0)} 
in the chiral limit in ( p5| ) and at temperatures near T* 
pion correlations are of the same size as f 2 (~ (a 2 ) in 
the 0(4) model) so that higher order corrections become 
equally important. 

We should bear in mind that the definition of f v is 
subtle even in thermal equilibrium. In fact, it is more 
convenient to define it as the residue of the axial-axial 
thermal spectral function < T c [A°(x), A^y)] > 
where A'tix) is the axial current (a = 1,2,3), instead 
of using the PC AC theorem jL7|. In this way, one avoids 
dealing with the reduction formula and asymptotic states 
at finite temperature. In addition, it is important to bear 
in mind that, due to the loss of Lorenz covariance in the 
thermal bath, one can define two different pion decay 
constants and /* corresponding, respectively, to the 
space and time components of the axial current In 
fact, the chiral symmetry imposes relations between them 
and the in-medium pion dispersion law. If the chiral 
symmetry is exact (as it is in our case) pions remain 
massless but their velocity u„ can be less than the speed 
of light and their thermal width can be different from 
zero. The relation with /*'* is given by v n ~ Re/*/Re/£, 
while the thermal width is proportional to the imaginary 
parts of /* and /* jH. Nevertheless, to one loop in 
ChPT one has Im/» = Im/* = and /« = /* = / ff (T) in 
( |30| ) although there are corrections to the velocity beyond 
one loop p3fl . As we are going to see, in the present model 
we will get a small but nonzero difference — /* to one- 
loop, unlike equilibrium, which could be interpreted as 
a small nonequilibrium deviation for the pion velocity. 
The pion velocity plays also an important role in the 
hydrodynamics of the chiral phase transition |29] . 



A. Nonequilibrium pion decay functions 

All the above considerations for /„■ can be extended 
to nonequilibrium. We refer to our earlier work |^8| for 
further details. The axial current from ([!]) reads 

AtfS, t) = *^tr [r a (U^d.U - Ud^)] (31) 

As in equilibrium, there are two i ndep endent f%(t) and 
fn(t). Their definitions are given in p8| consistently with 
the Ward identities of chiral symmetry. To leading or- 
der in ChPT (tree level with the action ([j])) one has to 
consider only 0(jr) contributions when expanding the U 
fields in the axial current ( |3l| ) which yields, as it should, 
/»(*) = /<(*) = f(t) at tree level. 

To next-to-leading order (NLO) there are three differ- 
ent contributions to f n : the first one comes from the 
NLO corrections to the propagator. Such corrections are 
of two types: one loop diagrams from the S2 action and 
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FIG. 3. The different diagrams contributing to the pion 
decay functions to NLO in ChPT 

tree level ones from S4 (diagrams a) and b) in Figure || 
respectively). It is important to bear in mind that in the 
calculation of diagram a) one has to integrate over all 
the branches of the contour C in Figure [l] and the result 
is independent of ti. Second, there is another one- loop 
diagram involving the product of three pion fields at the 
same space-time point (diagram c) in Figure ||) from the 
next order in the expansion of the axial current. Finally, 
at the same order we have to take into account that the 
axial current i n (Blj) is itself modified by the fourth order 
lagrangian in (fl6|) as [1 + fi(t)]Ao and [1 + fz(t)]Aj with 
/i(t) and/ 2 (i) in®. 

Once the different contributions have been taken into 
account, the final result to 0(p 4 ) reads j2^] 

lf°(t)} 2 = f(t) [1 + 2/ 2 (t) - h(t)} ~ 2iG (t) (32) 
[/<(*)] 2 = f(t)[l + f 2 (t)}-2iG (t) (33) 

for t > 0, with Ga(t) the equal-time correlation function 
defined in (13). The results fl3^)-fl33|) reproduce the equi- 
librium result in ( |30| ) when the time derivatives of f(t) 
are switched off and Gq 9 is replaced by j25|). 

As in standard ChPT, the loop corrections to f n come 
directly from the equal time pion correlator, which is 
time-dependent now. In addition, there are 0(p 4 ) tree 
level corrections given by the terms proportional to /1 (t) 
and f2(t). In the chiral limit in equilibrium there are no 
C(p 4 ) tree level corrections because Gg ? is finite. How- 
ever, at nonequilibrium Go(t) is UV divergent as we have 
seen in section [II and the counterterms proportional to 



fifiit) are precisely those needed to arrive to a finite an- 
swer. The above result for f^(i) should be such that the 
total answer is finite and scale independent because / w is 
an observable. Thus, we should be able to absorb the UV 
infinities and the scale dependence (in DR) in the n ew 
low-energy constant Ln, as discussed in section II C. In 
addition, from (p2|)-(|33|) and ( [l9| ) we see that the differ- 
ence [f%(t)] 2 — [/^)premains finite (it depends only on 
L12) so that the same renormalization is valid for 
and fl(t), which is another consistency check. In fact, 
the above result yields f%(t) ^ fl(t) to one loop, unlike 
the equilibrium case. Therefore, the plasma expansion 
induces modifications in the pion velocity larger than in 
equilibrium. However, note that we are following the 
equilibrium arguments given in fl3f in order to relate a 
nonzero value for /* — /* with the in-medium pion veloc- 
ity. Thus, our conclusions in this respect must be taken 
with care. In section [V C we will come back to this point 



B. f%(t) and fi(t) in parametric resonance. 

Let us concentrate now in the parametric resonance ap- 
proach, with f(t) given in ( pp| ) and where we keep only 
the leading order in the amplitude of the oscillations q. 
The loop contribution is given by Go(t), which we have 
analysed in detail in section HI. It will grow exponen- 
tially in time due to explosive pion production, once the 
infinities have been suitably subtracted, whereas the tree 
level corrections fi,2(t) remain bounded in time. 

The first step is to show how the infinities cancel in 
the final answer for f n (t). For that purpose, we replace 
in d3^)-([33"|) the functions fi,2(t) in (J19|) to leading order 
in q. On the other hand, in section [Fn] we have regular- 
ized the equal-time two-point function in DR. Its diver- 
gent divergent part for d — > 4 is given in (|29| ) . According 
to our previous discussion, we should be able to absorb 
the divergent part in In. In fact, collecting the piece 
proportional to L\\ in both /^'*(t) (remember that they 
only differ in terms proportional to L12) plus the diver- 
gent contribution in (^) and using (|l7|) yields 



[U(t) dlv ] 2 = -qM 2 sin Mt [12L\^) 



1 + log 



16tt 2 



4fc 2 



M 2 



(34) 



where we have taken the c? — *• 4 limit. The above con- 
tribution is finite and scale independent (the explicit de- 
pendence with [i is compensated by that in L\ x {\x) as 
explained before) which is a very important consistency 
check of our approach. Notice that it is crucial that the 
divergent contribution in (|29|) h as exactly the same time 
dependence as f(t)/f(t) in (|l^). 

Therefore, collecting the various pieces above, we can 
write the final result for f%(t) and f^(t) in parametric 
resonance (to leading order in q and to one loop in ChPT) 
as 



sin Mt 



1 



■ sin Mt 



L12 - 6In (fj.) 



qM 2 f I(t) 
ti \32tt2 
1 



32tt 2 



1 + log 



4k 2 



M 2 



- A B (t)}- 



t(t) 



J 7T 



(35) 



fl{t) = f°(t) (L 12 - -L 12 ) (36) 

for * > 0, with U(T) in @, I{t) in © and 

1 f°° , s sva.Mt- M: sin2kt , 

A s(t) = / dkkn B (k) ^- (37) 



8tt 2 



2vr 2 



dkk 1 



k 2 - M 2 /4 



iG (t,t,k)- — [l + 2n B {k)} 



and give some numerical estimates. 



(38) 
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Note that the only unstable (exponentially growing) 
contribution in (|3^) is given by A unst (t) in (p8|). The 
rest is bounded in time. Therefore, according to our pre- 
vious discussion, we can estimate the back-reaction time 
as A unst {tBR) — li- Thus, we have an approximate idea 
of the time scale during which our one- loop approach can 
be trusted. From that time onwards, the back-reaction 
corrections to f(t) coming from the coupling to the pion 
fields cannot be ignored. 



C. Numerical results 

Since the result is independent of the scale, we will fix 
fi = 1 GeV and use the numerical values of in and L12 
dete rmin ed phenomenologically in |34j and given in sec- 
tion |ll C[ . We have also taken ko = M (see our previous 
comments) . 

There are still three parameters we have to fix corre- 
sponding to the initial conditions: the initial amplitude of 
the oscillations q, the initial frequency M and the initial 
temperature Ti. According to our previous comments, 
M should be around the value of the a mass in the 0(4) 
model, although in the present approach it is not nec- 
essary to assume the existence of a a particle. We have 
considered three different cases: M — 0.1, 0.6 and 1 GeV. 
The second value is the one more often used in the litera- 
ture and also the closest one to the recent determinations 
39 1 . As for the initial temperature, we have fixed 
50 MeV. The standard approach 



of m 

for definiteness 7$ 
is to assume an initially supercooled (zero temperature) 
state when the a field starts its rolling down from the 
top of the potential. In this process the initial poten- 
tial energy is converted into thermal energy of the pion 
gas (reheating) |Q. Hence, it is reasonable to assume in 
our case a nonzero but small initial temperature. Even 
though our model does not have a direct interpretation in 
terms of an effective potential, it corresponds in the lan- 
guage of the 0(4) model to an initial condition where the 
a has already reached the bottom of the potential and it 
is oscillating around it. Nonetheless we must point out 
that the results depend very weakly on Ti. The reason is 
that in the right hand side of ( p5|) the dependence with 
Ti enters through the integrals (|37|)-(|38|) which are dom- 
inated by contributions near k ~ M and are therefore 
strongly damped by the Bose-Einstein distribution ns 
for temperatures Ti <C M. We have checked that tak- 
ing different values for Ti in the range 10-100 MeV the 
curves showed below remain almost unchanged. Finally, 
we have considered q — 0.1 and q — 0.2 to illustrate the 
dependence of our results with the initial amplitude. The 
initial values we are considering here are similar to those 
used in the literature |^, , M . 

The results for A unst (t)/f% are plotted in Figure ||. 
The values of the estimated tsn are also given in that 
figure. One clearly observes that increasing either q 
or M makes the unstable modes grow faster and over- 
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FIG. 4. The solid line is A unst (t)/f%. The dashed line is 
the tree level contribution l-q/2 sin Mt. T = 50 MeV for all 
cases. The values of the back-reaction time are approximately 
given by: (a) t BR — 300 fm/c, (b) t BR — 28 fm/c, (c) t BR — 
11 fm/c and (d) t B R — 10 fm/c 



come earlier the tree level value. In fact, the upper 
envelope of the long-time oscillations is proportional to 
(qM 2 /(4nUf)exp(qMt/2) since, from fl5§) and @), 
the dominant exponential contribution to the two point 
function at long times i s ex p(2/iz), the maximum of the 
Floquct exponent \l in (AE) being /1 ~ q/2 at the center 
of the unstable band, which has width q. 

The time it takes for the pion correlations to overcome 
the tree level has been estimated in the Q(4) model in 
this regime and it lies between 5-10 fm/c P,p6|. As ex- 
plained before, this is the same time scale as that when 
dissipation makes pion production stop. In the four cases 
we have considered, it is clear that our approach remains 
valid for the time relevant to pion production. For def- 
initeness, we will restrict from now on to the choice of 
parameters (d) in Figure ^ giving tsn — 10 fm/c, which 
is of the order of the plasma lifetime. In turn, note 
that for cases (c) and (d) we have gM 2 /A 2 ~ 0.07 and 



gM 2 /A 2 ~ 0.05 respectively, so that our ChPT approach 
to 0(q) is perfectly valid. 

In Figure |^ we have depicted the total result for 
f%{t)/ fir(Ti) in ( |35| ) for the same cases as in Figure [|. In 
this curve, the upper limit in time corresponds roughly 
to the onset of the back reaction. Its effect is to make the 
amplitude of the oscillations grow, whereas for t < tsR 
the amplitude remains approximately constant. As a 
matter of fact, one could wonder whether the oscillations 
of the unstable part could cancel those of the tree level 
in (|o|) so that fn{t) remains roughly constant at long 
times, which would be interpreted as a dissipation effect. 
Clearly, this is not the case as it can be seen in Figures |] 
and |f| The oscillations coming A unst (t) have indeed the 
same frequency as the tree level ones but their phase is 
different as it can be seen in Figure f|. In fact, in that 
figure one observes that for long times the phase of A unst 
is shifted almost 7r/2 with respect to the tree level and 
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therefore the contribution with cos Mt dominates. [] In 
other words, the effect of the pion correlations itself is 
not enough to make the system equilibrate. This is con- 
sistent with the analysis in J27],|37j where it is shown that 
such dissipation appears only when the back-reaction is 
self-consistently included, which we have not done, as ex- 
plained before. This issue will be confirmed by the anal- 
ysis in the next section, where it is shown that the pion 
number grows exponentially in time even when the loop 
corrections are included. If the back-reaction was con- 
sidered, the pion number should reach a maximum value 
and then stop growing [p7[ . Our f(t) is an external force 
whose shape is not changed during the time evolution, 
which is consistent only below tsR- Furthermore, not 
only the phase of the long-time oscillations is different, 
but, as it can be seen from Figure ^|, the central value 
also decreases with time (increases in Figure Q). This 
is due to the constant term proportional to exp(qMt/2) 
commented before. We interpret this effect as a reheating 
of the system [Q . The fact that decreases with tem- 
perature gives support to this idea. Also, assuming that 
when the system reaches the back-reaction time scale we 
can use approximately the equilibrium expressions then, 
according to ( pp| ) , an estimate of the final temperature is 
given by 



where the bar denotes time average: 

F(*j) = r r ' dt'F{t') 
t f Jo 



(39) 



(40) 



For the above estimate of Tf it is not important 
whether we choose ft or ft. Thus, for the parameters in 
case (d) in Figures |] and || we obtain, for Mt f ~ 25 — 30, 
Tf ~ 125 — 140 MeV which is not far from experimental 
determinations of the freeze-out temperature EJ|. Re- 
call that this value is almost independent of the initial 
temperature Ti and it is therefore compatible with a su- 
percooled initial state. Although our estimate is based 
on assumptions about the final state, we are using f n and 
not the pion distribution function. The reason is that 
always remains close to its equilibrium value, according 



^Working out the expressions for the unstable band given 
in Appendix |X| it can be seen that there are three differ- 
ent terms proportional to the leading long-time exponentials 
exp(gAfi/2) in Go(t,t,k). The first one is proportional to 
sin Mt and changes sign, to leading order in q, under a reflec- 
tion with respect to the band center. The other two are, re- 
spectively, time-independent and proportional to cos Mt and 
they are symmetric under such reflection. Therefore, by in- 
tegration in k, the sin Mt term is suppressed by at least 
one power of q with respect to the other two. The time- 
independent term increases the central value of the oscilla- 
tions, as seen in Figure W. 
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for the same cases as in figure 



to our previous discussion, unlike the distribution func- 
tion which is not thermal, as it will be shown in the next 
section. 

Next we will come back to the issue of the pion velocity. 
To the order we are considering, we have, from ( |35| ) and 



qM 2 
2Li2 sin Mt 



(41) 



It is unclear whether one can extrapolate from the 
analysis in Ej| and thus identify the above with the 
pion velocity u x . In fact, note that L\i < and then 
tin > 1 whenever sin Mi < 0. This is a similar problem 
as trying to identify m (t) with a pion mass which, as 
we have emphasized, is not correct since the chiral sym- 
metry is exactly preserved and the pions remain mass- 
less. The problem, as it is discussed below, is trying to 
define a time-dependent dispersion law, which is mean- 
ingless unless additional restrictions are imposed, such 
as adiabaticity. However, and following our previous ar- 
gument, by the final time where the equilibrium expres- 
sions are meant to be approximately valid, an estimate 
of the maximum variation of would be given by tak- 
ing the time average of the above quantity. Thus we get 
Av™ ax ~ 0.003 for Mt f = 30, for the same parameters 
as before. Note that this correction is even smaller than 
what is expected from ChPT beyond one loop, namely 
Avn ~ 0.13 for T = 150 MeV El. 



V. THE PION NUMBER 

A. The nonequilibrium particle number. Definitions 

One should bear in mind that the concept of particle 
number out of equilibrium is rather subtle. The nature 
of the problem is well illustrated once more in curved 
space-time QFT Q. The particle number depends on 
the reference frame and thus the initial vacuum state may 
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contain particles during its subsequent time evolution. In 
other words, the state which is regarded as the vacuum 
at time t is different from that at t = 0. However, the 
particle number can be given a physical meaning in some 
particular cases. For instance, if the metric is conformally 
Minkowskian (as the RW spatially flat metric we are con- 
sidering here) or it is Minkowskian on the space-time 
boundary. Another interesting regime is when the ex- 
pansion rate is small compared to the typical frequencies 
involved, which is the so-called adiabatic limit. Similarly, 
in nonequilibrium field theory one has to specify the state 
with respect to which particles are defined. One possi- 
bility is to choose the initial Fock space, which is sim- 
pler if the initial state is the equilibrium one (the analog 
of the Minkowski limit). In that case, one considers the 
time evolution of the initial creation and annihilation op- 
erators with the time-dependent density matrix p7| , ^7| . 
Physically, this corresponds to the number of initial par- 
ticles. In our case these are massless pions. Since our 
model preserves the chiral symmetry at all times, pions 
remain massless for t > 0. Given the difficulties related 
to the definition of the nonequilibrium dispersion law (see 
comments below), we will restrict here to the number of 
massless pions. Another possibility, when dealing with 
time-dependent mass terms like the one in (|) is to de- 
fine the adiabatic number, where the dispersion law is 



assumed to be simply uj 2 (t) = k 2 + m 2 (t) . The particle 
number is then defined in terms of the creation and an- 
nihilation operators that diagonalize the instantaneous 
Hamiltonian via a Bogoliubov transformation p3|| . We 
remark that the adiabatic limit is consistent for a slowly 
varying m 2 (t) and in fact it can be defined only for real 
Uk (t) . In our case this approximation would be valid only 
for k 2 > qM 2 /2. For small q this would capture anyway 
the resonance band and hence the parametric amplifica- 
tion. However, as we will see below, it is not even clear 
whether the dispersion law can be assumed to be adia- 
batic for all times. Both definitions of particle number 
coincide at t = 0. 

The nature of our approach makes it more suitable to 
define the particle number in terms of correlation func- 
tions, rather than in the canonical formalism. To illus- 
trate the way we will proceed, let us consider first a free 
scalar field <p of mass m in thermal equilibrium, so that 
n(k) is just given by the Bose-Einstein distribution func- 
tion ( |l2| ) . The time- independent Hamiltonian of the sys- 
tem is 



H = I d 3 x 



4> z + {Vcj)) 2 + m 2 4> 2 



(42) 



The thermal averaged energy of the system per unit 
volume is related to the particle number n{k) simply as 



d 3 k 



;(k) 



n(k) + - 



V J (2tt) 3 
where k = \k\ and u: 2 {k) = k 2 + m 2 is the free dispersion 



law. Thus, from p^ ) an d the definitions of the two-point 
functions in section II B . 



n(k) 



2uj(k) 



d d 



G>(t u t 2 ,k) 



ti=t 2 



dt\ dt 2 
oj 2 (k)G^(tx=t 2 ,k)] 
^— [-G>(0,k)+co 2 (k)G>(0,k) 



where Go is the free propagator, the dot denotes 
time derivative and we have used time translation in- 
variance (thermal equilibrium) so that Gq (t\,t2,k) — 
Gq (ti — t2,k), meaning that the particle number is time- 
independent in equilibrium. Note also that Gq above 
actually stands for G Q >11 (x,x') — ((f)(x)<f)(x')) since we 
are taking all time arguments in the C\ branch in Figure 
[| Using [d 2 + k 2 + m 2 ] G^ (t, k) = (the equation of 
motion) we have: 



i(k)=tu J (k)G>(Q,k)-l- = 



1 



1 



where we have used the solution for the free propagator in 
equilibrium which can be read off from jic| ) with hi(t, k) 
in (g), i.e, 

zG >e9 (i - f, k) = -L { [l + n B {k)] e- ik ^ 

+ n fl (fc)e ifc < t -*'>} (43) 

Our next step will be to extend the above definitions to 
the nonequilibrium case. Thus, following the same steps, 
we will define the nonequilibrium particle number n(k, t) 
through: 



W)) 
v 



= N 



d 3 k 
(2tt) 3 



to(k,t) 



n(k, t) 



(44) 



where oj(k,t) is (formally) the dispersion law, (E(t)) is 
the total energy of the system and N is the number of 
particle flavours. In our case N = 3, the number of 
different pions. In fact, we should consider a more general 
definition involving a sum over all internal indices but we 
will see below that to the order we are considering, all 
our expressions remain diagonal in isospin space. Thus, 
the pion density for a given pion type is 



V 



d^k „ , 
yPT n (M) 



(27T)° 



(45) 



where we have kept the space-time dimension d in or- 
der to regulate the UV behaviour (see below). Several 
remarks are in order here: As explained before, both 
the particle number and dispersion law may depend on 
time. In particular, loops can introduce corrections to 
the tree level dispersion law. Furthermore, whenever the 
time dependence appears through the interaction with 
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an external source (as it is the case here) the energy E(t) 
is not conserved (it is time dependent) as it happens in 
curved space-time when the back-reaction effects on the 
metric are ignored. Note that E(t) is the contribution 
to the energy from the pions only and the oscillations of 
f{t) transfer energy to the pions making the pion number 
n(t) grow with time. 

In fact, and following once more the curved space-time 
analogy we have previously discussed, we will calculate 
the expectation value of the total energy through: 



(E(t)) 



-0<2o (*>*)> 



(46) 



where T^ v is the energy-momentum tensor defined in 
terms of the lagrangian as p5[| 



(47) 



In Appendix [b| we have reviewed some useful results re- 
garding the calculation of the classical energy-momentum 
tensor for the case of interest here (spatially flat RW met- 
ric). We remark that, by construction, T Mt/ is symmetric 
and classically conserved, i.e, T" = where ; denotes 
the covariant derivative. However, this does not imply 
necessarily that the energy (Era) is time-independent. The 
reason is that P M = J d 3 x^/—q T.,° is not a covariant four- 
vector j}5| . Note also that in (|l6|) we are assuming that 
the expectation value of has perfect fluid form, which 
will be the case here (see below). 

Therefore, we will proceed by computing the energy- 
momentum tensor to a given order and then calculate the 
number of particles through ( f44| ) and (^) . This means 
to deal with expectation values of products of fields at 
the same space-time point and therefore divergent. We 
will follow the approach of point-splitting the fields so 
that the results are written in terms of Green functions. 
It is important to stress that, as it will be seen below, 
as long as we use dimensional regularization the final 
expressions for the particle number are automatically fi- 
nite, without any need for extra renormalizations of the 
energy-momentum tensor p0| , p5| , ^6[ . 

The simplest example of such point-splitting is pro- 
vided by our previous derivation of the equilibrium par- 
ticle number for a scalar theory, where we have replaced 
for instance: 



< cj> 2 (x, t) >= 



lira 



d d 



(x',t')->(g,t) dt' dt 



iG > (x-x',t-t') (48) 



We will proceed in the same way for the noncquilibrium 
case. However, as we will see in the next section, some- 
times we will have to deal with a field structure in the 
classical Too which is not symmetric under field exchange, 
like 4>(f>. In these cases we will symmetrize first the clas- 
sical expression and then point-split the fields J30| ] , i.e, 

(Ai(x)A 2 (x) . . . A n (x)) 



= i lim [(A 1 (x 1 )A 2 (x 2 ) . . . A n (x n )) 
+ (A 2 (x 2 )A 1 (xi)...) + ...] 



(49) 



where the Xj are space-time points, the Aj is a short- 
hand notation to denote either the field or an arbitrary 
number of its derivatives and the dots denote all possi- 
ble permutations. Note that we are dealing only with 
boson fields, which are symmetric under field exchange 
at different points. As it will become clear below, sym- 
metrizing the fields in this way yields consistent results 
for the particle number. 

Another problem in connection with the point-splitting 
is the Tc-ordering of the fields. For instance, we will find 
four-field contributions to the energy-momentum tensor 
and we need to relate them with T^-ordered four-point 
Green functions, so that we can use Wick's theorem to 
write the result only in terms of two-point functions. We 
did not have this problem in our previous derivation of 
the equilibrium particle number, since we were only deal- 
ing with products of two fields. Thus, we need to specify 
how the time arguments of the fields are ordered when 
taking the t, — > t limit. We will use the following pre- 
scription p5L|l] 

• First we will symmetrize over all possible ways of 
ordering the classical fields, as it is shown in (|49|). 

• Next, for a given ordering of the classical fields with 
an arbitrary number of field derivatives, we will re- 
place: 



= lim d Xl d X2 



(())(xi)(j>(x 2 )(f>(x 3 ) . . 
{Tcj){xi)(j){x 2 )(j){x 3 ) . 



(50) 



where lim*^^ means to take the Xj — ► x limit 
keeping the time arguments ordered from left to 
right, i.e, t\>t 2 >... (all tj are in the C\ branch of 
the contour in Figure [j] so that the contour ordering 
Tc becomes the ordinary time-ordering T). Note 
that the field derivatives, when present, arc pulled 
out of the T-ordering HQ. 

Note that in the second step we are choosing a particu- 
lar way to take the limit. If such limit exists, the answer 
should be the same regardless of the order. In this re- 
spect it is important to bear in mind that expectation 
values of products of fields (like G > (x, x') and its deriva- 
tives) have always a well-defined equal-time limit, unlike 
T-ordered products where one has to be careful with ill- 
defined expressions such as £(0), S'(Q) and so on, when 
taking time derivatives. For instance, suppose that we 
wanted to use the above prescription with (4>4>) to relate 
it with G(xi,x 2 ) = (T(j)(xi)4>(x 2 )}. Then, 
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m)m) = ^m)m+4>{t)kt)) 



— lim (dt 

2 t'-t v 



d' t )G> 



(t,f) = -\im(d t +d , t )G(t,t') 



(51) 



where for simplicity we have omitted spatial arguments, 
which do not play any role here, and lim* means taking 
the limit so that t > t' and then the last step in the above 
equation holds. Had we taken the limit keeping t' > t, 
the answer would have been the same by continuity of 
G > (t,t') and its derivatives. In fact, note that this is 
equivalent to take t and t' in C2 in Figure ^ which gives 
the same answer for G > as taking them in C\ in the t' — ► t 
limit |33). Of course, we could have written the result 
directly in terms of G > without specifying the order of 
the time arguments and the answer is the same either 
way. In other words, Wick's theorem is trivial for two- 
point functions. However, for four-point functions the 
above continuity arguments apply as well and hence this 
prescription will allow us to use Wick's theorem in the 
standard way for Tc-ordered products B. 



B. The number of pions in parametric resonance 

1. Leading order 

Let us start with the lowest order in ChPT. As ex- 
plained before, to lowest order it is enough to consider 
the lagrangian in ([l^) with the RW metric. Furthermore 
to leading order we only need the 0(tt 2 ) terms in that 
lagrangian. That is, the "free" lagrangian given in (||) for 
the 7r(x, i) fields. The reason why we can neglect 0(tt 4 ) 
terms to leading order when calculating (E) can be un- 
derstood in terms of Feynman diagrams. What we are 
doing is starting from diagrams with a given number of 
vertices of different types and closing them in all possible 
ways. For instance, at tree level there is only one point 
and hence one vertex. Thus, the contribution from the 
0(7T 2n ) is a n-loop closed diagram and hence according 
to Weinberg power counting theorem |]l5| it contributes 
as 0(p 2n ) in the chiral power counting. 

Therefore, using the form ( |l5| ) for the S2 action in the 
parametrization of the Tt a (x, t) fields, we have: 



Jc * 



Hence, the energy-momentum tensor ( |47| ) to lowest or- 
der reads simply 

T$ 2) = S M 7f°ft,7f° - lg^d a n a d a n a 



m pion fields. The above result is the standard kinetic 
term in curved space-time. Using the equations of motion 
to this order (g^d^Tr 11 ).^ = it is straightforward to 
(2 2) 

check that T^J is covariantly conserved. Besides, from 

(2 2) 

the particular form of T^J above it is not difficult to 
see that its expectation value has perfect fluid (diagonal) 
form, that is, (Ty) = for i =/= j and (Tjo) = 0, just from 
spatial translation invariance. 

Now, let us consider the above for the RW metric. The 
total energy defined in (Eg) reads to this order 



« 2 W r_a l2 >(*) a^a 



(E^ (t)) = a 2 (t) *)] 2 + [V* a (x, t)f 

= i J <fx([Ti a ] 2 + [V7T a ] 



where we have used (B4) and the superscript (n, m) in 
T^ v means a contribution coming from the S n action with 



2— T 7T"7T' 

a(t) 



where in the last line we have written the result for the 
Tr a (x,t) fields. The first line above simply states that 
the energy is conformally equivalent to the Minkowski 
(equilibrium) result in the tt parametrization. Recall that 
the spatially flat RW metric is related to the Minkowski 
one by a conformal transformation. 

Note also that the energy density is not obtained 
just by replacing m 2 —> m 2 (t) for a free scalar the- 
ory, as one could have expected from the lagrangian 
(||). That would be equivalent to work in the adiabatic 
limit and it would have been the answer defining the 
energy as H(t) = J cPxTL(x,t) with the Hamiltonian 
density H — ^(dC/dir) — C. However, these two defi- 
nitions arc not equivalent in the presence of the external 
force a(t). In fact it can be checked that adding to the 
lagrangian a total derivative y/—gC — * y/—g (C + Afy 
with A M a contravariant vector (the action has to re- 
main a scalar) which is a functional of the field and 
its derivatives, does not change neither the equations 
of motion nor the energy-momentum tensor defined as 
@. However, H may change under such transforma- 
tion. For instance, considering the expression for the 
action S^tt] before integrating by parts to get (0), i.e, 
5 2 [tt] = (l/2)/(a p 7r^7r + (d/a) 2 7r 2 -2(d/a)7r7r), one 
gets H = (1/2) (-k 2 + (Vtt) 2 - (d/a) 2 7r 2 ). Therefore in 
that case the "mass" term would be proportional to 
— (a/a) 2 rather than —a/a. In this case A = (a/a)ir 2 . 
As we have said before, none of these expressions for the 
energy density is time-independent. 

Nevertheless, there is a way to check the consistency 
of the above result which in fact gives us a hint of how to 
include the back-reaction effects. If we let a(t) (or f(t)) 
be a classical field (independent of x for simplicity) in the 
lagrangian, then, by considering also its equation of mo- 
tion, it can be checked that (E^ 2 ^(t)) — 0. Moreover, 
H(t) coincides with ( |52"|) when the extra term coming 
from Ti — > H + a(dC)(da) is considered. From this point 
of view, in the present approach we are just consider- 
ing that the field a does not receive quantum corrections 
and its equation of motion is dominated by the kinetic 
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term, which depends only on g^ v but not on the pion 
fields. That term is nothing but the counterpart of the 
Einstein-Hilbert action, yielding the Einstein equation 
for when no matter fields are present ||3q] . 

At this point let us recall that we are considering only 
the leading order in q consistently with our power count- 
ing. Therefore we can neglect the d 2 term in (|52|). Now, 
let us apply our point-splitting prescription to the re- 
maining terms in (52). We will write the result in terms 
of hi(t) in ©. Let us consider first the term propor- 
tional to d. According to our previous discussion, 



(ir a (x)n a (x)} = 3i lim 



r d 3 k 1 


' d d~ 


/ (2tt) 3 2 


dt + ~dt' 



d 3 k 
(2^)3 



2n B (k)]^-\h 1 (t,k\ 2 
dt 



Si 



Gq (t, t' , k) 



G (t) (53) 



where we have used fllCQ and (|14|). The above result de- 
serves some comments: first, we realize that if we had not 
taken the symmetric limit, the answer would have been 
complex in general, whereas it is manifestly real when the 
two contributions are added together. This is a consis- 
tency check of our point-splitting prescription. Second, 
since hx(t, k) = iexp(— ikt)/y2k + 0(q) (the leading or- 
der is the equilibrium solution) then (TT a -k a ) — O(q). In 
other words, for the time scales we are considering (below 
t B R) Go{t)/f% — 0{q). Therefore and since d = 0{q) as 
well, to leading order we just have 



{E^ 2 \t)) 



V 



d 3 k 
h!(t,k)\ 2 



{[1 



2n B (k)} 

fc 2 M*,fc)l 2 



}+0{q 2 ) (54) 



Hence, according to (fMJ), we find that the number of 
massless pions of a given type is given to leading order 
by 



n(k, t) 



^ [Mi^P + fc^iM)! 2 ] [l + 2n B (k)]-l 

(55) 

This expression coincides with the one given in [ ^7|j24| 
in the canonical operator formalism. Note that to this 
order the answer for the adiabatic number of pions would 
have been the same, since m 2 (t) = 0{q). Note also that 
our initial conditions (JsJ) imply that n(k, 0) = n B (k) as 
it should. 

So far we have considered the general result in ChPT, 
for arbitrary f(t). Let us now particularize to the para- 
metric resonance case in (p0[). First, in order to analyse 
the UV behaviour of the pion density in , we need the 
solution h% to the Mathieu equation in the stable band 
(large k). Using the approximate solutions to leading 
order in q (see Appendix |X|) one has: 



qM 2 sin Mt - M- sin 2 kt , 0s 
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FIG. 6. n(k, t) for M = 600 MeV, T; = 50 MeV, q = 0.2. 



\h 1 (t,k)\ 2 = - 1 



qM 2 sin Mt - f|sin2/tf 
~~4 k 2 - M 2 /A 



-0{q 2 ) 



for k > fco (see section III). Thus, we see that the possible 
logarithmic divergence in d — > 4 (i.e, that coming from 
the term which is not proportional to n B {k) in j55|)) can- 
cels exactly to leading order in q and therefore the total 
number is finite. This justifies our approach of point- 
splitting the fields as long as we remain within the DR 
scheme. 

The pion number will grow exponentially in time due 
to the contribution of the unstable band. By the same ar- 
guments as those used for the two-point function, n(fc, t) 
will grow in time typically with exp(qMt/2). In Figure 
6 we have plotted n(k, t) for the choice of parameters (d) 
in Figure [|. We observe that the pion number grows to 
order one within the resonance band and before the back- 
reaction time. We have also plotted the time average of 
the pion distribution function n(k, t) in Figur e ^ taking 
Mt, 



30 (see our comments in section IV C| ). 



We observe the typical peak of parametric resonance 
amplification at k ~ 300 MeV, which would be seen in the 
final pion spectra pj| although single pion distributions 
might be not enough to disentangle from states which are 
not DCC-like and one needs to consider higher order pion 
correlation functions Note that the final shape is 

very different from a thermal Bose-Einstein distribution 
function, although both diverge for k — > 0, which re- 
flects Bosc-Einstein condensation for massless particles. 
The physical pion number density in momentum space 
is k 2 n(k,t) 1 whose time average we have also plotted in 
Figure [7] (dashed line). 

Finally, the pion density (|^|) as a function of time is 
showed in Figure |[ Our results agree numerically with 
the predictions of the 0(4) model in the spinodal regime 
[^,^2|. As we have previously commented, the particle 
number grows in time without stop. One may wonder 
whether the one-loop NLO effects may change this pic- 
ture, but this is not the case, as we are going to see in 
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FIG. 7. The solid line is the averaged pion number n(k,tf) 
for t f = 30/M and the dashed line is (k 2 /M 2 )n(k,t f ). Here, 
T, = 50 MeV, q = 0.2 and M = 600 MeV. 




t [fm/c] 

FIG. 8. The pion density for T t 
M = 600 MeV. 



50 MeV, q = 0.2 and 



the next section. The origin of this behaviour, as empha- 
sized before, is that we have not taken into account the 
back-reaction properly |2?j . As long as we remain below 
the back-reaction time, our results for the pion number 
can be trusted. 



2. Next to leading order 

There are three types of NLO contributions to the pion 
number: 

1. The one-loop corrections to the two-point fun ction, 
which we have already discussed in section IV A 
(diagrams a) and b) in Figure ||) . 



2. The contribution of the 0(p 4 ) lagrangian (|l6| ) to 
the energy-momentum tensor. According to our 
previous discussion, only two-field terms contribute 
to NLO, namely (E^ 4 ^) in our notation. 

3. Four-field terms in the 0(p 2 ) lagrangian, i.e, 
(£( 2 > 4 )). 

Before proceeding, we should insist that, as com- 
mented before, we are considering the number of mass- 
less pions, ignoring possible one-loop modifications of the 



dispersion law. We must stress that the difficulties as- 
sociated to the definition of the time-dependent disper- 
sion law are similar to those related to the nonequilib- 
rium particle number. Physically it makes sense to de- 
fine it asymptotically at long-times but it is not clear 
whether one can actually define an instantaneous disper- 
sion law, unless one follows the adiabatic approximation. 
Nonetheless, we believe that such corrections would not 
change qualitatively the time evolution of the pion num- 
ber. In section pAJ we have already commented on the 
in-medium dispersion law and its relationship with the 
pion decay constants fp . Remember that in equilibrium 
the pion dispersion law remains unchanged to one-loop 
in ChPT. While a detailed nonequilibrium extension of 
the results in is out of the scope of this work, qual- 
itatively we expect a similar one-loop behaviour in our 
case with perhaps a change in the pion velocity of order 
!/*(*)//* (t) | as (|l|) shows. That contribution is pro- 
portional to L±2 and bounded in time. In any case, if we 
write u> = k+T,(k, t) with S(fc, t) = O(q) a bounded func- 
tion in time, the contribution to the pion number ( |44| ) 
would be n — ► n — (S(fc, t) / k)riB{k). This is a bounded 
correction that does not change the relevant features of 
pion production we are analyzing here. 

Let us then start with the NLO corrections of type 
[j]. According to our power counting, this correction af- 
fects only the terms in E^ 2 - 2 \ The NLO correction to 
the pion two-point function has been calculated in |2§| l 
as explained before. Let us write it as G > (t,t' , k) = 
Gq (t,t',k) + A(t,t',k). Then, according to our argu- 
ments in the previous section, the contribution of this 
correction to the pion number is given to 0(q) by: 



n(k, t) — > n(k, t)+ 



2k 



d__d 

dtdt' 



7 A(t,t',k) 



+ k 2 A(t,t,k) 



Now, from the results in |2^| (eqs (15)-(17) in that 
paper *) we obtain: 

i[d t d t ,A(t,t',k)\ t=t , + k 2 A(t,t, k)] 

- -ih(t)k 2 G>{t,t) - l -h{t) [d t + d t ,]G>{t,t')\ t=t , 

+ [2A 2 {t)-ih{t)] d t d t ,G>{t,t>)\ 



duAxi 



+ k 



, Gi 



duA? 



(u,k) {[(d t G>(t,u)) 2 - {d t G>(u,t)) : 
*(t,u)) 2 -(G>(u,t)f}} 
(u) { [(d u d t G>(t,u)) 2 - (d u d t G>(u,t)Y 



* There are two misprints in eq.(16) in J28|. The term 
6f(t)/f(t) multiplying G (i) should read 4f(t)/f(t) and the 
term —2Go(t) should read — Go(t). None of them affects the 
results here, since they are included in the Ai(t,k) function 
in (IBl). 
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(d u G>(t,u)) 2 -(d u G>(u,t)) 2 



(56) 



Here, the fc-dependence of Gq (t,i? , k) has been sup- 
pressed for simplicity, is the function in ( |l9| ) ap- 
pearing in the C(p 4 ) lagrangian and we have used that 
d t [G$(t,t') - G£(t,t')\ = -1. The explicit form of the 
functions Ai(f, k) and Aa(t) is given in p8[ . What is 
important for our purposes here is that Aj.,2 = 
within the range of validity of our approximation, i.e, 
while G (t)//| = O(g). In fact, A 2 (t) = G Q (t)/f 2 (t), 
while the precise form of Ai is unimportant here (see 
below). Hence, to 0{q) it is enough to replace in ( |5q ) 
the leading order term in q for Gg" , which is nothing but 
the equilibrium solution (p3|) . [j] This simplifies consider- 
ably the above expression. In fact, the dependence with 
Ai(t, k) disappears and (p6[) reduces to: 



A(t,t',k) 



dtdt 
-k[l + 2n B (k)} 



+ k 2 A(t,t,k) 



</ 7T 



+ 0(g 2 ) (57) 



Therefore, the correction to the pion number from this 
part is, neglecting 0(q 2 ), 



n(k,t) ^ n(k,t) - ~[l + 2n B (k)} 



ZiW+i-p- 

Jtt 



(58) 



where / x (t) = -12 (2L U + L 12 ) /(*)/./* + 0(q 2 ). 

Note that in our one-loop calculation of f v in section 



IV 



we have shown that the combination —YlL\\fjf — 
iGg(t) is finite. This is the combination appearing in 
(|32|)-(]33|) with fi(t) and / 2 (t) in ([l9"l). However we have 
here —24Luf/f + iGo(t), which diverges. The only pos- 
sible way out is then that the remaining NLO corrections 
(types ^] and ^|) combine with this one in such a way that 
the answer for the pion number is finite. 

Let us then consider type || corrections, i.e, those com- 
ing from the energy-momentum tensor to fourth order 
in derivatives. By the same argument as before, only 
two-field terms contribute to this order. Hence, from the 
lagrangian in ( [l6| ) we can concentrate only in the L\\ and 
L12 terms. The energy-momentum tensor coming from 
that lagrangian is calculated in Appendix |b|. The final 
expression for T^t' 2 ^ is displayed in (B7) for arbitrary 
metric and the contribution to the energy for the RW 
metric is given in (B9). Note that we could have used 
the equations of mot ion to second order to simplify some 
of the terms in ( |B9| ), writing for instance tt in terms of 
Air and tt 2 and so on. It is clear that to 0(q) we only 



^This holds also in the unstable band without expanding in 
q in the leading exponentials exp fj,Mt where [i = O(q) is the 
Floquet characteristic exponent (see Appendix |a|). 



need to consider the terms proportional to g\ , 55 and 57 
in that expression. We have 



= — lim 



( ^ k 2 


' d d~ 


1 w 


dt + dF 



G>(t,t',k) (59) 



and 



(7r a (x,t)n a (x,t)) 



= — lim 

2 t'^t 



r d 3 k d d 


' d d~ 


) {2nfdtdt' 


db + dt 1 



G>(t,t',k) (60) 



where we have retained only the LO propagator Gq 
since the two terms above are multiplied by g§ (t) and 
gj{t) which are already 0(q). Now, since Gg (i, t',k) = 
G> eq {t - t',k) + 0{q) with GT q {t- t',k) in @ it is 
clear that the contributions (p9h-(pQ) are both O(q) and 
therefore they do not contribute to the order we are con- 
sidering here. 

Hence, the only correction of this kind is the one given 
by the term proportional to gi(t) in (| B9 | ) . T his term is 
of the form already analyzed in section~ |V B 1 . From our 
results there, we have that the correction to the particle 
number from (E 1 - 4 - 2 ^) is then given to O(q) by 



n(k, t) -> n(k, t) + [l + 2n B (k)] 



2gi(f) 

f3 

J 71 



(61) 



where gx{t) - -3{2L n + L 12 )f(t) + 0{q 2 ). 

Finally, we shall consider the contribution of four fields 
from the 0(p 2 ) lagrangian in (fLq) to the energy momen- 
tum tensor. Expanding to 0(W^) in that lagrangian we 
find : 



T, 



(2,4) 



00 



3 

■' a,b=l 



T^a^b [^a^b + (V7f a ) (Vfffe)] (62) 



Now we should write the energy ( [46|) to this order for 
the 7r° = TT a f(t)/f fields, from the above expression. 
However, according to our power counting, since this con- 
tribution to the energy is 0(p 2 ) with respect to the tree 
level, it must be at least 0(q). This will be confirmed 
by our subsequent calculation. Therefore, when chang- 
ing from the tt fields to the tt fields we can simply ignore 
the terms proportional to f(t). Hence, to this order it is 
enough to replace simply tt by tt and f(t) by / in ( [62] ) 
and therefore 



E {2A) j d 3 xl Ka(x,t)TT b (x,t)[TT a {x,t)TT b {x,t) 



0,6=1 

+ (Vir a (x,t)){V7T b (x,t))]} 



(63) 



We have to calculate now the expectation value of the 
above quantity. Note that this is the first time where we 
find the problem of the T-ordering discussed in section 
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V A , apart from the issue of the symmetrization of the 



classical fields. We will follow the prescription explained 
in that section for the different field structures appearing 
in (|63|). Let us consider first the terms ('K a 'Ki ) d 0i 'K a d a 'K b ) 
with a ^ b and where the Lorenz index a is fixed, i.e, is 
not summed over. We have 

(7r Q (x)ir b (x)d a ir a {x)d a ir b (x)) 

= — lim [d^d^ i (TTT a (x 1 )7r b (x 2 )7r a (x 3 )Tr b (x 4 )) + . . .} 

4! %j — >x 

where the dots stand for all the permutations of the four 
fields in the expectation value. Now, using Wick's theo- 
rem, 

{Tit a (xi ) 7r b (x 2 ) ir a (x 3 )^b {xa ) ) 

= -Go(x u x 3 )G (x 2 ,x i )[l + 0(p 2 /A 2 x )] for a^b 

which is the only order we have to retain according to 
our power counting (remember that this contribution to 
the energy is already NLO and therefore it is enough to 
keep the LO propagator) and we have taken into account 
that the lowest order action in (j|) is diagonal in isospin 
space, so that to lowest order the T-product of fields with 
different indices factorizes. Note also that the result to 
this order is independent of the a, b indices as long as 
a b. Collecting the different permutations we find 



(ir a (x)Tr b (x)d a ir a (x)d a Tr b (x)) 
~ {dl+di)G>{x lX ') 



-i 2 



(a^b) 



It is clear, following the same arguments as before, that 

the above quantity is 0(q 2 ) since (d x + Gq(x, x') — 

O(q). Thus, we only need to consider the terms with 
a = b in ( |63| ) to 0(q). For a given n field and a fixed, 

(ir 2 (x)ld Q ir(x)] 2 ) = l lim 



6 xj 



+ d Xl d X3 + 



(T7r(xi)7r(:E2M:E3Mx4))- 

Thus, Wick's theorem gives now 

(T7r(a;i)7r(a;2)7r(a;3)7r(x4)) = - [G (xi, x 2 )G (x 3 , x±) 
+ G (xi,x 3 )G {x2, xi) + Go(x 1 ,x 4 )G (x 2 ,x 3 )} 

Now, we take the Fourier transform in the spatial com- 
ponents of the above expression. Our prescription for the 
lim* Xj ^ X IS equivalent to repl ace a ll the Gq above by 



Gq (see our comments in section V A ). Thus, taking into 



account once more that Gq (t,t' , k) = Gq q (t — t',k) + 
O(q), G>(t,t,k) = Gp(t,t,k) and G (t)Jf 2 = 0(q), we 
find for the combination appearing in (p3|), 



iG {t) 



n 2 (x,t){l7r(x,t)] 2 + {Vn(x,t)} 2 }) 
d A - x k 



Therefore, from (|63|) we find that the contribution of 
the corrections of type [| to the particle number is given 
by 



n(k, t) -> n(k, *) + [! + 2n B (k)] 



iG (t) 



(65) 



(2tt) 



-k[l + 2n B {k)]+0(q 2 ) 



(64) 



Finally, collecting the contributions to the pion number 
to NLO, namely, (p8|), ( |6l| ) and ( |6q ) we find an interesting 
result: the total NLO correction to 0(q) vanishes. Note 
the completely different origin of these three contribu- 
tions and remember that each of them was UV divergent 
so that there were only two alternatives: either they can- 
cel or they appear in the combination (32)-([33|). Thus, 
the only NLO correction to the pion number is due to 
the change in the dispersion law. Note that, in turn, we 
have shown the absence of NLO corrections to the pion 
number for equilibrium ChPT. This is indeed consistent 
since we know that in equilibrium in the chiral limit the 
pion dispersion law is unchanged to one loop in ChPT. 
Therefore, the pion distribution function has to be the 
Bose-Einstein one for massless particles since the system 
remains in thermal equilibrium. 

Thus, the numerical results showed in Figures 6 and| 
remain valid to NLO. As commented before, we did not 
expect that the NLO corrections make the pion number 
stop growing since we have not included the back reaction 
and the energy E(t) to this order is still not conserved. 
Our results are valid below the back reaction time and 
should account for all the relevant pion production. 



VI. CONCLUSIONS AND OUTLOOK 

Chiral perturbation Theory can be used to describe 
nonequilibrium phenomena. In particular, in this work 
we have showed that pion production can be accommo- 
dated in ChPT in the parametric resonance regime. The 
physical situation where this analysis is meant to be use- 
ful is the late time expansion of the plasma formed after 
a Relativistic Heavy Ion Collision. Pion production is 
important in the context of hadronization and produc- 
tion of Disoriented Chiral Condensates during the chiral 
phase transition. 

In the present approach, we have considered the Non- 
linear Sigma Model in the chiral limit, where the pion 
decay constant is time-dependent. This is a nonequi- 
librium effective model with a well-defined pcrturbative 
expansion and power counting near equilibrium. Besides, 
using the analogy of this model with curved space-time 
QFT, we have been able to construct the fourth order la- 
grangian and implement renormalization in a consistent 
fashion. The parametric resonance regime corresponds 
to take fn(t) oscillating around its equilibrium position. 
To lowest order f^t) corresponds to the vacuum expec- 
tation value of the a field in the 0(4) model. Thus, the 
pion equation of motion to lowest order in the amplitude 
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oscillations becomes a Mathieu equation, which has res- 
onance bands in momentum space. The pion correlator 
grows exponentially in time, yielding explosive pion pro- 
duction. This approximation is consistent until the time 
when the back reaction effects due to the pion correla- 
tions become important. We have estimated this time 
scale for different choices of the initial values of the am- 
plitude, frequency and temperature. 

The two observables we have analyzed here are the 
pion decay constants and the pion number up to one 
loop in ChPT. Our main results are the following. For 
fn{t) the nonequilibrium corrections are basically of os- 
cillatory nature until the back reaction time. However, 
the central value tends to decrease, which can be inter- 
preted in terms of a reheating of the system. Besides, a 
small difference between /* and f* is induced, unlike the 
equilibrium case where it vanishes at one loop. Using the 
equilibrium result, we have estimated the final tempera- 
ture and the averaged value of — ]\ which is related 
to the in-medium pion velocity in equilibrium. 

As for the particle number, we have first introduced 
a suitable definition in terms of the energy-momentum 
tensor and Green functions. We have showed that in di- 
mensional regularization there is no need for extra renor- 
malizations and we can use a point-splitting prescription 
consistently. The number of initial particles at tree level 
coincides basically with the result in the 0(4) model. At 
one loop we have found that all the relevant contribu- 
tions cancel, to leading order in the oscillations ampli- 
tude. This result holds also in equilibrium, which is a 
particular case of this analysis. Thus, our prediction for 
the particle number is just the tree level result, which 
gives pion exponential growth in time. We have given 
numerical results both for the pion distribution function 
n(k, t) and for the pion density (n(t)) /V. The final distri- 
bution function which would be observed has the typical 
peak of parametric resonance at the center of the unsta- 
ble band (k ~ M/2). 

The reason why f n (t) and (n(t)) /V are not damped is 
because we have not taken into account the back-reaction 
effects which would change the original ansatz for the 
pion decay constant. When taken into account, those ef- 
fects should make the particle number stop growing and 
give energy conservation. Nevertheless, our approach is 
perfectly valid until the time where these dissipation ef- 
fects are important and therefore we believe we capture 
the essential behaviour concerning explosive pion produc- 
tion. 

We must stress than, apart from being a physically 
interesting case, the example analyzed here has allowed 
us to show explicitly the renormalization of our model 
to one loop, which is not trivial because of the presence 
of new nonequilibrium infinities. In fact, we believe that 
our methods could be useful for other nonequilibrium 
field theoretical models. 

There are many directions in which this work can be 
extended. Perhaps the most important would be to be 
able to include the above mentioned back-reaction effects, 



in order to understand dissipation properly. Other exten- 
sions include to consider nonzero physical pion masses 
and other relevant observables such as the correlation 
length or higher pion correlators, which are important to 
clarify the issue of DCC formation and to obtain predic- 
tions testable in RHIC. In addition, photon production 
can be studied by gauging the NLSM, including the Wess- 
Zumino-Witten term, responsible for the anomalous de- 
cay 7T° — > 77. Work along these lines is in progress. 
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APPENDIX A: SOLUTIONS OF THE MATHIEU 
EQUATION 

Here we will summarize the main results used in 
the text concerning the solutions of Mathieu differential 
equation. All these results can be found in ]47|j4^ |. 

According to Floquet's Theorem, there is always a so- 
lution of the Mathieu equation ( ]2~]] ) of the form 

F v {z) = e wz P{z) (Al) 

where P(z) is a periodic function with period tt and v is 
called the characteristic exponent, which depends on a 
and q and it plays a crucial role in our analysis, since it 
gives rise to exponentially growing solutions whenever it 
takes complex values. 

The values of a such that F v {z) is periodic in z are 
called the eigenvalues of the Mathieu equation. They 
correspond to integer values of v. They are denoted as 
a r (q) if v is a positive integer r, and b r (q) if v = —r. 
It can be shown [|^J that for a > one has a r > b r , v 
is complex in the bands b r < a < a r and real elsewhere. 
Moreover, for small q, one has b r —a r = 0(q r ). Therefore, 
in the narrow resonance regime we are considering here, 
we will take v complex for b\ < a < a\ and v real for 
< a < b\ and a > a\. The series expansion in q of the 
eigenvalues is given by 

ai = l + q + 0(q 2 ) ; b x = l-q + 0{q 2 ) (A2) 

If a ^ a r ,b r then F v (—z) is a solution linearly inde- 
pendent of F u (z). This is no longer true if v is an integer, 
although an independent solution can also be constructed 
in that case pgj ]. For v 2 7^ r 2 , it is customary to take as 
independent solutions: 

ce v {z) = \[F v (z)+F v (-z)\ 

se v {z) = ^[F v {z)-F v {-z)\ (A3) 
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For v = r, ce r (z) are called the eigenfunctions of the 
Mathieu equation and so on for v — —r and se r {z). They 
are 27r-periodic and their g-expansion for v 2 — 1 is given 
by[||: 

cei(z) = cosz — — cos(3z) ; se±(z) = sin z — — sin(3z) 
8 8 

(A4) 

Therefore, the solution h\{z,k) of ( pl| ) is given as a 
linear combination of ce„ and se v : 



h\{z, k) — A{k)ce v {z, k) + B(k)se v (z, k) 



(A5) 



where the coefficients A(k) and B{k) are such that the 
initial conditions dO) are satisfied, i.e, 



fi = -V(oi -a){a-h) 



C\ = \J a - bi+ ] 



S± = \/ ai — a — ji 2 



(A9) 



Note that fi — 0(q) and it reaches its maximum value 
at the band center. The dominant behaviour at long 
times is therefore given by the positive exponentials when 
(|A|) is replaced in @ and @. 



APPENDIX B: RESULTS IN CURVED 
SPACE-TIME 

1. General results 



A(k)ce v (-%/4, k) + B(k)se v (-ir/4, k) = 
A(k)ce v (-n/4, k) + B(k)se v (-Tr/4, k) = 



2k 
2k 



AI 



where the dot means d/dz. Remember that ce„ and se v 
depend on k through a{k). 

Even though the solutions to the Mathieu differential 
equation are numerically tabulated, we need their explicit 
form when dealing with renormalization. Such a explicit 
form of the solutions can be found for small q. Let us 
consider first the case when v 2 ^ r 2 and real, i.e, the sta- 
ble zone. Then, using the g-expansions given in p7| , p8[ , 
one has v = \fa + 0(q 2 ) and the solutions are given by: 

ce v (z) — cos(y/az) + q [cos(y / az) cos2z 
+ Vasin(\/az)sin2z] +0(q 2 ) 

se v (z) — sin(v / a2) + q \sm(-\/az) cos2z 
— \/a cos^^/az) sin2z] + 0(q 2 ) 



(A6) 



where 



1 



2(a- l) 2 -q 2 



(A7) 



We will collect here some of the results concerning 
curved space-time needed for our purposes. In this sec- 
tion we will consider an arbitrary metric g^ u and in the 
next one we will particularize for the spatially flat RW 
metric. Most of the definitions used here can be found in 
any textbook on the subject (we are following the nota- 
tion and conventions of fl35|). 

The covariant derivative of a contravariant vector 
V fl (x) satisfies 

V=gv.£ = 9 M (V=^ M ) => J d'x^VZ = (Bl) 

which is the generalized Gauss theorem and V^(x) is as- 
sumed to vanish at the space-time boundary. 

The Ricci tensor and scalar of curvature are defined 
respectively as R^ v 
Riemann tensor is: 



R\ v and R = g^R^ where the 



R%& — ds^p^ — d^Tpg + r^ A r /37 — r^r^ (B2) 

and the Christoffel symbols are given in terms of the 
metric as 



1 



9 Xa [dn9va + d v g^ a - d a g^ 



(B3) 



If a is far enough from the border points of the first 
band (placed at a — 1 ± q) we can simply take q ~ 
q/[2(a — 1)] in (A6). The value of a from which this 
simplification is valid can be estimated numerically, for 
a given q, by comparing to the numerical (tabulated) 
solutions and impos ing that the difference with the ap- 
proximate solution (A6) remains 0(q 2 ). 

Now, consider the un stab le band, i.e, b\ < a < a\. In 
this case, the solution (Al) reads, to leading order in q 



\C\ce\(z) + Sisei(z)] 



(A8) 



where ce\{z) and se\(z) are given asymptotically in (A4) 
and the real characteristic exponent \i and the Ci, Si 
coefficients are given by 



The classical energy-momentum tensor of the the- 
ory T^^x) is defined by performing a general coordi- 
nate transformation (infinitesimal) g^ v (x) — > g^(x) + 
Sg fiu (x), under which the action S = J y/—gC changes as 
S —>■ S + 8S where by definition 



SS =2 



^gT^Sg^ 



The energy-momentum tensor thus defined is symmet- 
ric and conserved (T^. v = 0) as long as the action is a 
Lorenz scalar. Under the above transformation the vari- 
ation of the metric determinant is given by 



(B4) 
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whereas the variation of the Ricci tensor yields 



SRuu = \g 



v;X 



P V '; P \X 



(B5) 



which is known as the Palatini identity. 

In the text we need the energy-momentum tensor de- 
fined in ( p7|) , to fourth order, i.e, when the lagrangian is 
given by (|16|) . Since we are interested only in two-point 
functions, we can ignore the contribution of C±(U, g) and 
concen trat e only in the Lu and L%2 pieces. Therefore, 
using (B4), 



T^(x) = -g, v (x)L^\x) + 2 S 4^M (B6) 



5g^{x) 



with 



C^ 2 \x) = - [L 11 R(x)g^(x) + L 12 R^(x)} F^(x) 
F pu = tr [d^(x)d u U(x)] = ^d^ a (x)d„Z a (x) + 0{^) 

in the parametrization of the tt(x) fields. 

To calculate the variation of the lagrangian in (B6) we 



use (B5) and integrate by parts taking into account (Bl) 
as well as the properties of the covariant derivative. We 
find after a straightforward but lengthy calculation: 

T$ 2) = -Lu [{2R^ - 9fiU R) g af) F a(i + 2RF„ V 

— 2g P vg al3 Faf3' S $ + g^ {F a p. <tl . v + Faftv-p) 

— L\2 [2 (R^Fva + R^Ffia) — g f _ l vR a(3 F a p 



g^F a0 ^F^ s +F0.^ + FP, 



V ;0 



(B7) 



The above energy-momentum tensor is symmetric, 
since g M „, R^ v and F^ v are symmetric. Note that the 
symmetry of F^ v is a consequence of U being unitary so 
that 



On the other hand, in Minkowski space-time, where 
R P v = and the covariant derivatives are ordinary par- 
tial derivatives, we recover the result given in |34f| . Fi- 
nally, we have verified explicitly that T^ 4 ' 2 ) = 0, using 
the equations of motion for £i 4 > 2 \ which to 0(tt 2 ) read 

{[LnRgi™ + L 12 R^} = (B8) 



2. Results for the RW conformal metric 

As explained in the text, the space-time metric we are 
interested in is the RW spatially flat metric in conformal 



time, where the scale factor is a(t) = /(£)//• The line 
element is ds 2 = a 2 (t)[dt 2 — dx 2 ] so that the elements of 
the metric are: 



ffoo (*) = a 2 (t), gij{t) = -Sija 2 (t), g i0 = 

and the metric determinant is g = detg — —a s (t). The 
nonvanishing Christoffcl symbols for this metric are: 



P o /.\ _ a (t) r k _ x k r o _ >- 
m[) ~a(t)' °*~a(t) 1 ' Lij ~^) dij ' 



o(t) 
a(t)' 



d(t) 



and the nonvanishing elements of the Ricci tensor are 
given by: 



Roo(t) = 3 
so that the scalar of curvature is 

a(t) 



"a(t) 


a 2 {tY 


, R^ — 


~a(t) 


, « 2 W1 


a(t) 


a 2 (t)_ 


a(t) 


a 2 (t)\ 



R(t) = 6 



a*(t) 



With the above ingredien ts we can calculate the 
energy-momentum tensor in (B7) for this metric. It is 

easy to check that (T^' 2 ') is diagonal for the RW met- 
ric above, i.e, (T- ( 4 ' 2) ) = and (T^' 2) ) = if i ^ j, as 

(2 2") 

w i^n,,,,,,. oi so w ith the lowest order (T^ ) (see sec- 



tion 


VB1 


tota 


ener 



VB 1 ) . Therefore, we only need Too to calculate the 



result here in terms of the ir fields and the f(t) function, 
retaining only two-field terms: 



d 3 xa 2 (t)T$ 2) 

! {gi(t) [*°] 2 + <? 2 (*) [W Q ] 2 + g 3 (t) [tt-] 



E^ 2 \t) 
4 

+ g A (t)-K a * a + g 5 (t)(X7Tr a ) ■ (Vtt q ) + g 6 (t)n a n a 
+ g 7 (t)rr a n a } 



(B9) 



where we have integrated by parts neglecting total spatial 
derivatives, the space-time dependence of the fields has 
been suppressed for simplicity and 



(?L n + 2L t2 ) 
-3(5L n +L 12 ) 



fit) 2 



gi{t) = -3(2L U + L 12 )f(t) - 3(5L n + L 12 ) 

92® 

g 3 (t) 

54 (t) 



fit? 
/(<) 



/(*) 
fit) 4 



,/(*) 



/(*) 



f 3 (t) 

(2L n +L 12 )f(t)+2(5L n 



L12 



fit? 
fit) 
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gs(t) = -(6L 11 + Lu)f(t) 
g 6 (t)=-3(2L n + L ia )^r- 
g 7 (t) =3(2Lu+Li 2 )/(t) 



(BIO) 
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